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Abstract 

One of the goals of the Large Hadron Collider (LHC) is to search for super symmetric particles 
predicted by the minimal Supergravity Model (mSUGRA). In previous studies the discovery 
potential for the next-to-lightest neutralino (X2) has been investigated in many areas of its 
parameter space. It has been pointed out that the leptonic decay of the X2 offers kinematic 
properties which can be well used for mass reconstructions of sparticles. However, in the param- 
eter space where taus dominate the leptonic decay, the X2 discovery is especially difficult since 
a full tau reconstruction will not be possible at LHC. In several recent publications it has been 
indicated that studies are necessary in that region. 

This study analyses the possibility of using kinematic endpoints for a sparticle mass recon- 
struction at tan/3 = 35, Aq = and \x > 0. A region of 25 points around the benchmark point I' 
with mo = 181 and mi/2 = 350 is chosen where the cascade decay X2 — ► r ± + f T — ► r ± + t t + Xi 
has a branching ratio of over 95%. In this cascade decay at least three particles, the Xi and 
two f T , escape the detector. Therefore, a precise mass reconstruction in the investigated region 
using the kinematic limits is not possible. However, this study shows that some invariant mass 
endpoints can be well estimated, using a mo and mi/2 independent correlation between the dis- 
tribution maximum and the kinematic limit. Furthermore, a linear fit at the tail of the invariant 
mass distribution gives in most of the cases a good estimate as well. In addition, a method is 
described how to measure two endpoints within the single rho with one quark invariant mass 
distribution. The influence of the mixing between b\ (pi) and the heavier 62 (£2) on endpoint 
measurements is also discussed. 

The simulations are performed with PYTHIA 6.220 in combination with ISASUGRA 7.69. 
The results are at Monte Carlo level, and give the minimal systematic and statistical uncer- 
tainties on measurements of 13 kinematic endpoints. For the mass reconstruction of sparticles 
with a 30 fb' 1 data sample the minimal uncertainties of 20.4 GeV (7.9%) for the X2, 20.5 GeV 
(14.8%) for the Xi, 21.2 GeV (14.0%) for the ff , 28.2 GeV (3.6%) for the d x and Si, 25.7 GeV 
(3.6%) for the 62 and 28.8 GeV (3.9%) for the £2 have been calculated. 



Chapter 1 

Introduction 



Today, the Standard Model [1] of electroweak and strong interactions describes nature at the 
smallest scales accessible in high energy physics. Both theories are based on the same funda- 
mental principle, the local gauge invariance. Through precision measurements at LEP and SLC 
the electroweak theory has been checked very thoroughly, and no obvious deviation or incon- 
sistency has been observed. The only particle of that theory that still has to be discovered is 
the Higgs boson. It is necessary for explaining the existence of massive fundamental particles in 
nature. Quantum Chromodynamic (QCD), describing the strong interactions within the Stan- 
dard Model, is tested very well at HERA. Only for gravitation, it has not been possible yet to 
construct a theory based on quantum mechanics which can be tested in experiments. 

However, there are a lot of open questions which have to be solved. The origin of mass, the 
reason for the huge matter antimatter asymmetry in the universe and the nature of cold dark 
matter are not known until now. The Standard Model has a large number of parameters which - 
instead of arising from a more fundamental theory - have to be measured. Furthermore, there are 
theoretical motivations for physics beyond the Standard Model. One is the hierarchy problem 
concerning the quadratically divergent fermion loop corrections to the Higgs mass. Additionally, 
- for reasons of concinnity - theorists want to unify the gauge couplings. J.C. Maxwell working 
out his famous equations as well as S. Glashow, S. Weinberg and A. Salam - all three responsible 
for the electroweak unification - have shown that a unification of theories can be very successful. 

In order to answer some of these questions, one possibility is to extend the Standard Model 
Poincare algebra to a supersymmetric graduated Lie algebra which leads to a larger particle 
spectrum. The model with the minimal particle content is the Minimal Supersymmetric Stan- 
dard Model (MSSM). It is an important prediction of this theory that every particle described 
by the Standard Model has a supersymmetric partner, a sparticle: Each fermion has a bosonic 
and each boson has a fermionic superpartner. However, none of the predicted sparticles has been 
discovered yet. If our universe really was supersymmetric directly at the Big Bang, Supersym- 
metry had to be broken at a very short time after it. Otherwise, every particle would have the 
same mass as its superpartner, which is definitely excluded by experiments. The MSSM solves 
the hierarchy problem, provided there are sparticles in the 1 — 10 TeV mass region. It makes 
a unification of the gauge couplings possible, and in some cases even predicts a candidate for 
cold dark matter. Furthermore, the MSSM gives strong constraints for the lightest Higgs mass 
which there has to be lower than 130 GeV [2]. 

The mechanism for breaking Supersymmetry is not specified in the general MSSM. The 
breaking through gravity is one possibility, and is given by the minimal Supergravity model 
(mSUGRA) [3]. This theory is very appealing because it reduces the 124-dimensional parameter 
space of the MSSM to five dimensions. Since in mSUGRA an ideal candidate for cold dark matter 
can be predicted, also experimentalists use this model. 

The start of the Large Hadron Collider (LHC) at CERN near Geneva is foreseen in the 
year 2007. In addition to searches for the Higgs boson, one of the main goals is to test the 



predictions of prospective theories. If the MSSM describes nature, it is the aim to produce 
enough supersymmetric particles such that the two detectors CMS and ATLAS are able to 
discover them. After the discovery, it is important to measure their properties since these 
properties are necessary for distinguishing between different ways of Supersymmetry breaking. 
It is thus essential to determine the best methods for sparticle searches within the different 
models before the start of the LHC. 

A short introduction into the principles of Supersymmetry is given in chapter [21 There the 
focus is set on the Minimal Supersymmetric Standard Model (MSSM), since it provides the 
simplest supersymmetric extension of the Standard Model being consistent with experiments. 
Furthermore, the mSUGRA model is introduced being one possible realisation of the MSSM 
and the framework in which this study has been carried out. The Monte Carlo generators used 
in this work are presented in chapter |3 In chapter @J three important experiments concerning 
this study are presented. Part of the mSUGRA parameter space has been excluded due to limits 
on sparticle masses with the Large Electron Positron collider (LEP) at CERN. The Wilkinson 
Microwave Anisotropy Probe (WMAP), a NASA satellite, provides strong constraints for the 
cold dark matter density in the universe. In mSUGRA the Xi is the lightest supersymmetric 
particle. Assuming the Xi to be the candidate for cold dark matter, the parameter space can be 
further substantially reduced. In the future, the Large Hadron Collider (LHC) will be an ideal 
machine to search for Supersymmetry. Chapter |S] gives an introduction to the basic kinematics 
used in this analysis. The formulae and configurations for the kinematic endpoints are explained. 
Chapter is the main chapter where the analysis methods and the results for X2 searches at 
high tan/3 are presented. Many surveys have been made in several areas of the mSUGRA 
parameter space which are accessible at the LHC. It has been shown that the leptonic decay of 
the next-to-lightest neutralino, 

X°2 -» f* + l T -» l± + l T + Xi , 

has a useful kinematic feature: The dilepton invariant mass spectrum has a sharp edge near 
the kinematic upper limit. This feature was first discussed in [4]. If the decay into electrons 
and muons has a sufficient branching ratio, the endpoint can be measured very precisely, and 
therefore can be used for sparticle mass reconstructions. In some regions of the parameter 
space, however, the decay into taus dominates. At the benchmark point I' with tan/3 = 35, 
Aq = 0, (jl > 0, mo = 181 and m\/2 = 350 [5] 96% of all X2 decays are into taus. Therefore, 
methods introduced for the dileptonic decay of the X2 are applied to the tau channel. Since 
the endpoint is always smeared out due to the presence of neutrinos in the final states, the 
existence of a correlation between the different endpoints and the distribution maxima has been 
investigated. For distributions where the ratio between both values shows a linear dependence, 
the measurement precision can be improved. The single rho with the associated quark channel 
offers two endpoints which can be well estimated in most of the cases without knowing whether 
the corresponding tau comes from the slepton or the X2. Additionally, it turned out that by 
using more endpoints than sparticle masses the determination of the sparticle masses can be 
improved significantly. The outlook in chapter |H1 summarises problems and ideas which have not 
been treated in this work. 



Chapter 2 

Supersymmetry 



The Standard Model successfully passed all particle physics experiments in the last 30 years. 
Especially the LEP and SLC experiments have tested its predictions with high precision. Indeed, 
there are experiments like the measurements of the anomalous magnetic moment of the muon 
in Brookhaven [6] and the WMAP measuring the cold dark matter density in our universe (see 
section I4.2JI which may give a hint that physics beyond the Standard Model exists. But they 
do not exclude it. Thus, the motivation for introducing Supersymmetry is based on theoretical 
aspects. 

In nature it is often possible to explain its principles and make predictions using symmetries. 
Since the discovery of the Pauli Principle, particles are divided in two groups, the fermions whose 
spin is an odd multiple of tj and the bosons with integer spin. A supersymmetric operation on 
the representation of a particle belonging to one group transforms it to a member of the other 
group. Thus, the spin changes but the mass and the couplings to other particles remain the 
same. The electron with spin i and a mass of 0.511 MeV would have a bosonic partner with 
spin and the same mass in a supersymmetric world. However, no supersymmetric partner has 
been discovered until now. It is therefore certain that - if Supersymmetry exists at a high energy 
scale - it is broken at the electroweak scale which is at the moment accessible for experiments. 

2.1 MSSM 

Particles and their interactions are described by Quantum Field Theories. The fundamental 
principles of nature are described through particles and their interactions whereas also the 
interactions are treated as particles in such theories. A quantum field is a complex n-dimensional 
function ip(x) representing one or more particles. tp(x) has the property that |t/>(£)| 2 gives the 
probability to find that particle at the space-time point x. Since only |^(x)| 2 is measurable, the 
phase of ip(x) is a free parameter. It is natural to assume that at any point in the universe the 
phase can be chosen differently without changing the particle's physical state: 

i>(x) -► e ieix) 4>(x). (2.1) 

This is the basis of all gauge theories like the Standard Model and its minimal supersymmetric ex- 
tension. The direct consequence is that gauge bosons are responsible for the interactions between 
fermions. In the Standard Model the electroweak interactions result from a local SU(2)l x U(1)y 
symmetry and the strong interactions from a local SU(3) symmetry. The interactions within 
the MSSM are based on the same symmetries. 

The MSSM is minimal in the sense that it is the simplest possible supersymmetric extension 
of the Standard Model with a minimal particle content. The MSSM has an additional second 
Higgs doublet, describes more than twice as many particles than the Standard Model and solves 
the following problems: 



1. One of the main problems of the Standard Model is the Hierarchy Problem. The radiative 
corrections to the Higgs mass including the fermion loops at one-loop-level has a unphysi- 
cal quadratic divergence. Indeed, this divergence can be canceled by a mass counter term 
which, however, has to be fine tuned at each order in perturbation theory. In the MSSM, 
two additional scalar particles S, provided they are in the region of 1 — 10 TeV, are re- 
sponsible for the cancellation of that divergence. The remaining divergence is then only 

i| - roj) 



logarithmic, multiplied with (m% — ml) being zero in the supersymmetric limit. 



2. Within the Standard Model, it is not possible to achieve a unification of gauge couplings 
- at any scale - which reduces the number of open parameters, and is therefore favoured 
by many theorists. The MSSM provides the possibility of gauge coupling unification at a 
scale of about 10 16 GeV. 

3. In some models within the MSSM there exists a candidate for cold dark matter. Experi- 
ments measuring the rotation velocities of galaxies have shown that it is very likely that 
there is more than just visible matter in our universe. One component of the so called dark 
matter seems to be cold, meaning that it consists of particles moving at a non-relativistic 
speed. Indeed, the Standard Model axion can be responsible for cold dark matter. But, it 
has not been discovered until now, and it is not clear whether it can explain the measured 
cold dark matter density sufficiently [7]. If mSUGRA is the model describing the Super- 
symmetry breaking in nature and if the lightest supersymmetric particle is a neutralino, 
it would be an excellent candidate for cold dark matter. 

4. In the MSSM the mechanism of electroweak symmetry breaking occurs radiatively at the 
level of perturbation theory, without the need for any new strong interaction. Thus, in 
contrast to the SM the MSSM gives an explanation for the Higgs mechanism. 

However, the MSSM gives neither a hint to the source of CP violation nor is it able to make 
predictions for any particle mass. Furthermore, the MSSM still is a theory without quantum 
gravity. String theories may solve some of these problems, although a phenomenologically viable 
string theory has yet to be constructed [8]. 

There are two different kinds of superfields in the Minimal Supersymmetric Standard Model 
at the GUT-scale and above: 

1. The Chiral Scalar Superfields consisting of a complex scalar field, S, and a 2-component 
Majorana fermion field, £. The SM quarks and leptons and the MSSM higgsinos at 
the TeV-scale are connected with £ and its partners, the squarks, sleptons and Higgses, 
with S. 

2. The Massless Vector Superfields consisting of a gauge field F^ v and a 2-component Majo- 
rana fermion field, Xa- The field F^ v is connected with the Standard Model gauge bosons 
and Xa with its Superpartners, the gauginos, at the TeV-scale. 

These fields respect the same SU{2>) x SU(2)l x U(1)y gauge symmetries as do the Standard 
Model fields. In order to confirm or confute the predictions of the MSSM it is first necessary to 
find the new particles meaning to find the mass peaks of the mass eigenstates (Column Mass in 
Table l2~Tj) . In the MSSM every supersymmetric particle has the same gauge couplings than its 
Standard Model partner. Thus, the second step would be to measure the interaction couplings 
and the spin of the new particles. 



MSSM Particles 


Extended Standard Model Super symmetric Partners 


Interaction 


Mass 


Spin 


Interaction 


Mass 


Spin 


e ± u* r ± 
^L/R^L/Ri ' L/R 


e ± ,/i ± ,r ± 


1/2 
1/2 


e* u* f* 


e* u* f ± 

c l/2'^l/2' 'l/2 






U L/R' C L/R' t L/R 
^L/_R' S L/_R' °L/R 


u,c,t 
d, s, b 


1/2 
1/2 


U L/R-> C L/R-> t L/R 
"■L/R^L/R^L/R 


^1/25^1/2,^1/2 






B°,W° 
H ui H d 


h,H,A°,A ± 


1 



B°,W° 

H"0 rrO 

H u^d 


Xl,X2,X3,X4 


1/2 
1/2 


w ± 


W ± 


1 
1 


W ± 


Xi,xf 


1/2 
1/2 


9 


9 


1 


9 


5 


1/2 



Table 2.1: This is the particle spectrum of the MSSM. The extended Standard Model has 
instead of only one an additional second Higgs doublet. The labels "L" and "R" denote the 
left- and right-handed electroweak eigenstates. The mass eigenstates Xi/ 2 are defined as a?i/ 2 = 
x l/r ' cos 9 ± Xrjl ■ sin9. Since 9 is close to zero for electrons, muons, up, down, charmed and 
strange quarks, the mass eigenstates are mostly denoted with "L/R" instead of "1/2", in the 
literature. 



2.2 Superpotential and R-parity 

In the SM, the effects due to the masses of quarks and leptons are parametrised by the Yukawa 
potential. The Yukawa couplings in the MSSM are specified by means of a function called 
superpotential. The general superpotential of the MSSM contains terms where the baryon 
and the lepton numbers are violated. These terms have significant experimental restrictions 
especially from measurements of the proton decay and the double beta decay. There are two 
ways to make the theory describing the data. First, the involved parameters can be tuned such 
that they are in accordance with experimental results. Alternatively, by introducing a new 
symmetry this problem is solved. This symmetry and its corresponding transformation is well 
explained in [8] and [9]. It leads to a multiplicative quantum number, the R-parity, which by 
convention is 1 for all Standard Model particles and —1 for their superpartners. It is given by 



R 



_^3(B-L)- 



-2,s 



(2.2) 



for a particle with baryon number B, lepton number L and spin s. This symmetry has an 
important influence on the MSSM phenomenology since it requires that supersymmetric particles 
always are produced in pairs and that no supersymmetric particle can decay into solely Standard 
Model particles. In particular, the lightest supersymmetric particle is stable. 

2.3 mSUGRA 



In mSUGRA it is assumed that the three interactions are unified at the Grand Unification 
Theory (GUT) scale (10 16 GeV) and that Supersymmetry is broken due to gravity. The beauty 
of mSUGRA lies in the fact that this model needs only four more parameters and the sign of 
a fifth additional parameter with respect to the known free parameters in the Standard Model. 
If mSUGRA describes the breaking of the Supersymmetry correctly, all scalar particles like 
sfermions and Higgs bosons have a common mass mo at the GUT scale. Also the gaugino 



masses M±, M2 and Ms, corresponding to the U(l), SU(2) and SU(3) gauge symmetries, unify 
to a common mass mjft at the GUT scale which can be seen in figure |2~T1 
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Figure 2.1: Evolution of the particle masses in mSUGRA through the Renormalisation Group 
equations (RNG) as a function of the energy scale. The mass spectrum is reduced to three 
distinct values at the GUT scale, whereas the particle spectrum at the experimentally accessible 
scale is split. In mSUGRA the gluino is always the heaviest gaugino. 



All trilinear coupling parameters A^ have the same value Aq at Mqut-, and the ratio of 
the vacuum expectation values of the two Higgs doublets is given by tan/3. The magnitude 
of the Higgsino mass parameter \i is given by the four mSUGRA parameters but its sign is 
the fifth open parameter. In mSUGRA with conserved R-parity (see section 12.2(1 the lightest 
supersymmetric particle (LSP) has to be stable. In some parts of the mSUGRA parameter space 
the 7\ is the LSP. However, a charged stable supersymmetric particle is cosmologically excluded. 
In other parts, the lightest neutralino is the LSP. This neutral weakly interacting particle is an 
ideal candidate for cold dark matter. If the latter is a concentration of lightest neutralinos, the 
WMAP data (see section fO|) . giving constraints on the cold dark matter density in the universe, 
can be used to restrict the mSUGRA parameter space. Very often, this restriction is shown in 
the rao — mi/2 plane for different values of tan/3, Aq and sign(//). 

One of the challenges of the LHC is to decide whether mSUGRA is the right theory for the 
description of nature. If this is the case, the next step will be to determine the values of the five 
parameters. 



Chapter 3 

Monte Carlo generators 



At LHC, protons - which are made of quarks and gluons - will be accelerated to 7 TeV, and 
will collide at different interaction regions. These collisions will produce leptons together with 
quarks and gluons, both fragmenting into jets, in the final state. In most of the cases it is not 
possible to calculate analytically the Quantum Field theory based production cross sections and 
4-momenta of particles produced in such collisions. Therefore, Monte Carlo generators, generally 
used for statistical models in science, are a useful tool for simulations in particle physics. In 
order to estimate the sensitivity of analysis methods, which will be utilised after having collected 
the first real data, simulations of proton-proton collisions at the LHC are necessary. Here, the 
first step is to generate quark-quark, quark-gluon and gluon-gluon collisions with a Monte Carlo 
generator. The particle distribution functions which describe how the quarks and gluons are 
distributed within the protons are measured at HERA and by other experiments. The choice 
of the special Monte Carlo generator depends on the physics model for which sensitivity studies 
are made. The output at this level are the precise 4-momenta of all particles in such an event. 
Then detector effects and effects due to event reconstruction have to be included. 

3.1 ISASUGRA 

In this study ISASUGRA 7.69 is used which is part of ISAJET 7.69 [10]. It calculates nu- 
merically the masses, the branching ratios and the weak-scale mass parameters within the 
mSUGRA model. The particle spectrum calculated by ISASUGRA can be used as an input 
for PYTHIA [11]. The cross sections are evaluated by ISAJET or PYTHIA. ISASUGRA takes 
the following values as input: The choice of the mSUGRA parameters and the top mass. The 
top mass is set to 175 GeV in this study which is the same value used by most of the recent 
publications of mSUGRA studies. Presently, the updated value is 174.3 ± 5.1 GeV [12]. Fur- 
thermore, Aq is set to zero, fi is positive and the value for tan/3 is 35. The values for mo and 
mi/2 and the motivation for the parameter choice are given in chapter |BJ 

3.2 PYTHIA 

PYTHIA is a program for the Monte Carlo generation of high energy physics events arising 
from collisions between elementary particles such as e + , e~, p and p in various combinations. 
It is based on theories and models for a number of physics processes, including hard and soft 
interactions, parton distributions, initial and final state parton showers, multiple interactions, 
fragmentation and decay. The version used in this study is the corrected (See appendix I A. 1|) 
version of PYTHIA 6.220. In PYTHIA it is possible to include new particles, cross sections, 
processes, and also to specify the decays of particles. In order to save CPU time, the taus in this 
study have been forced to decay solely into charged rhos. The production cross section calculated 
by PYTHIA did not change, with respect to the cross section with all tau channels included, 



due to that adjustment (See appendix I A. 4(1 . Therefore, in order to normalise all distributions 
to 30 fb' 1 , the PYTHIA 6.220 branching ratio of 24.9 % for r ± -► p ± v T has to be included. 
In this study an event is defined as one produced X2 decaying into two r. Thus, the number of 
produced events has been multiplied by the squared branching ratio. 



Chapter 4 



Experiments 



4.1 Large Electron Positron Collider 

The operation of the Large Electron Positron (LEP) collider at CERN started in August 1989. 
With a circumference of 27 km, LEP was the largest accelerator in the world. It stopped 
in November 2000, but the analysis of data is still going on. LEP accelerated electrons and 
positrons - in opposite directions in a vacuum pipe inside a retaining ring of magnets - which 
collided head-on at different interaction points. The four LEP experiments were called ALEPH, 
DELPHI, L3 and OPAL. 

The MSSM predicts a light Higgs with M h o < 130 GeV, and with LEP a lower constraint 
of Myfi > 114 GeV was obtained. In some parts of the mSUGRA parameter space, especially 
for low values of ttiq and mi/2, LEP should have produced enough super symmetric particles 
necessary for a discovery. However, there was no evidence for Supersymmetry at LEP. Thus, 
a value for tan/3 lower than 2.4 has been excluded [13]. The LEP Susy Working Group has 
obtained preliminary results which are shown in figure I4~T1 and H~2l These results are based on 



1. a combination of E C m = 183 - 208 GeV ALEPH, DELPHI, L3 and OPAL searches for 
pair production of e in the e — > eXi channel, 



2. a combination of E CM = 183 - 208 GeV ALEPH, DELPHI , L3 and OPAL searches for 
pair production of f in the f — ► rX\ channel, 



3. a combination of the E C m = 205 - 208 GeV ALEPH, L3 and OPAL searches for pair 
production of Xx , 



4. the Eqm < 202 GeV ALEPH searches for pair production of heavy stable charged particles, 

5. the DELPHI and ALEPH searches for neutralinos decaying into f and r, 



6. a combination of E C m < 209 GeV ALEPH, DELPHI , L3 and OPAL searches for the 
lightest scalar neutral Higgs in the ee — > hZ reaction. 



ADLO Preliminary - |U>0, M top =175 GeV/c , A Q =0 




200 400 600 800 1000 



200 400 600 800 1000 




200 400 600 800 1000 

m n (GeV/c 2 ) 



200 400 600 800 1000 



m (GeV/c ) 



Figure 4.1: The excluded regions in the mo-m^-plane through measurements with LEP for 
[i > and Aq = 0. Yellow (1): no mSUGRA solution: no EWSB or tachyonic particles; Light 
blue (2): regions inconsistent with the measurement of the electroweak parameters at LEP1; 
Green (3): regions excluded by chargino searches; Red (4): regions excluded by selectron or stau 
standard searches; Dark Blue (5): regions excluded by the search for hZ; Brown (6): regions 
excluded by the neutralino stau cascade searches. Magenta (7): regions excluded by the search 
for heavy stable charged particles applied to staus [14]. 
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ADLO Preliminary - |U<0, M top =175 GeV/c , A Q =0 
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Figure 4.2: The excluded regions in the mo-m^-plane through measurements with LEP for 
[i < and Aq = 0. Yellow (1): no mSUGRA solution: no EWSB or tachyonic particles; Light 
blue (2): regions inconsistent with the measurement of the electroweak parameters at LEP1; 
Green (3): regions excluded by chargino searches; Red (4): regions excluded by selectron or stau 
standard searches; Dark Blue (5): regions excluded by the search for hZ; Brown (6): regions 
excluded by the neutralino stau cascade searches. Magenta (7): regions excluded by the search 
for heavy stable charged particles applied to staus [14]. 
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4.2 Wilkinson Microwave Anisotropy Probe 




Figure 4.3: A view of the WMAP satellite in the universe [15]. 



The cosmic microwave background (CMB) radiation is the radiant heat left over from the Big 
Bang. It was first observed in 1965 by Arno Penzias and Robert Wilson at the Bell Telephone 
Laboratories in Murray Hill, New Jersey [16]. The properties of the radiation give information 
about the physical conditions in the early universe, and can therefore be used for a better 
understanding of the formation of cosmic structures such as galaxies and dark matter. As a 
consequence, assuming that cold dark matter is made of neutralinos, the CMB radiation data 
leads to restrictions on the allowed mSUGRA parameters. Thus, this data is important for 
making constraints on the mSUGRA parameter space. 

If the Xi were abundantly produced directly after the Big Bang they would have had mostly 
annihilated to the relic density which has been measured with WMAP. In neutralino relic density 
calculations for mSUGRA there is a 'bulk' region of relatively low values of TOj/2 and too, where 

neutralinos annihilate dominantly through t-channel slepton exchange (XiXi — ► I I). For a large 
neutralino mass, the relic density is generally large. However, neutralino-stau coannihilation 
(Xif — > tV) can reduce the relic density to satisfy the cosmologically favored relic density 



Mj-o <C M^o ) has a highly fine-tuned relic density 



bounds. The region where this is possible (M, 

since a small variation in too would lead to a large change in SlLSph 2 . Furthermore, there is a 
region where rapid annihilation through s-channel Higgs boson resonances (XiXi — ► H — > bb, tt) 
can reduce the neutralino relic density to the measured region. Finally, in the focus point region 
at large too - where the Xi is Higgsino like - the relic density is achieved mainly through the 
XiXi -» W+W~ and the XiX? -» ZZ channels. 
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Figure 4.4: The shaded strips display the regions of the (mi/2,mo) plane that are compatible 
with 0.094 < QLSph 2 < 0.129 at the 2-c level where Q,LSph 2 is the normalised density. These 
areas are in the bulk (B',C',G',I',L'), coannihilation tail (A',D',H',J') and rapid-annihilation 
funnel (K',M') regions. Left plot: fi > and tan/3 = 5 (A'), 10 (B',C), 20 (G'JT), 35 (I',J') 
and 50 (L',M'). Right plot: fi < and tan/? = 10 (D') and 35 (K'j. The parts of these WMAP 
lines for /i > compatible with the g^-2 measurements at BNL [17] at the 2-cr level have lighter 
(pink) shading [5]. 

The Wilkinson Microwave Anisotropy Probe (WMAP), shown in figure 14.31 is a NASA 
satellite, and presently the best experiment to measure the CMB radiation and its spatial fluc- 
tuations. Launched in 2001, WMAP is now able to resolve small temperature fluctuations of the 
order of only millionths of a degree. Instead of measuring the absolute temperature values of 
the CMB radiation WMAP measures the temperature differences between two different points 
in the sky, taking advantage of the higher accuracy of differential measurements. Figure 14.41 
shows that the region in the two-dimensional mi/2-Trio plane consistent with experimental data 
is reduced to nearly one dimensional strips due to the WMAP measurements. Therefore, the 
WMAP experiment gives a hint where we have to focus the searches for mSUGRA. 



4.3 Large Hadron Collider 

This study concerns physical processes that might happen at the Large Hadron Collider (LHC). 
The construction of the LHC has been approved by the CERN's Council in December 1994. The 
LHC, shown in figure 14*31 is the successor of the Large Electron Positron collider (LEP), and will 
be installed in the 27 kilometre circumference LEP tunnel being about 100 metre underground. 
The LHC will operate at a center-of-mass energy \fs of 14 TeV being the highest ever reached 
in the world. The design luminosity is 2.3 • 10 34 cm~ 2 s~ 1 . Furthermore, lead-ion collisions are 



27 



cm 



The accelerator 



planned with a total y/s of 1148 TeV and a luminosity of 1.0 • 10 
consists of two interleaved synchrotron rings, whose main elements are 1232 superconducting 
NbTi dipole and 392 quadrupole magnets operating in superfluid helium at a temperature of 
1.9 K. The dipole magnets guide the particles along the ring while the quadrupole magnets 
are focusing the particle bunches. The nominal field for LHC magnets to handle 7 TeV beams 
is 8.34 Tesla, and the goal is to achieve 9.0 Tesla [19]. Two 450 GeV beams coming from 
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Figure 4.5: A schematic view of the LHC at CERN near Geneva [18]. 

the Super Proton Synchrotron (SPS) will be injected into the LHC, and then accelerated with 
superconducting cavities, 8 cavities per ring, operating at 2 MV and 400 MHz. 

The experimental program of the LHC will include the dedicated heavy-ion detector ALICE 
(A Large Ion Collider Experiment), the specialised B-physics spectrometer LHCb and the general 
purpose detectors CMS (Compact Muon Solenoid) and ATLAS (A Toroidal LHC Apparatus). 
The physics goals of CMS (Figure I4.6JI and ATLAS include searches for the Higgs boson, which 
is predicted by the Standard Model, as well as for supersymmetric particles which appear in 
Supersymmetry theories like mSUGRA or the Gauge Mediated Susy Breaking (GMSB) scenario. 
If the latter particles exist they will be produced mainly through squark or gluino production 
at LHC. These processes can lead to long decay chains which allow to give constraints for all 
involved sparticle masses. 
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Figure 4.6: A schematic view of the CMS detector. The pixel detector is closest to the interaction 
vertex. The second innermost detector is the silicon tracker followed by the crystal electromag- 
netic calorimeter and the hadronic calorimeter. Outside of the superconducting solenoid which 
generates a magnetic field of 4 Tesla there are the muon chambers. The total weight is 12500 T 
which is nearly twice as much as the weight of the Tour Eiffel. 
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Chapter 5 

Kinematics for X2 searches 



At the Large Hadron Collider the colliding parton's center of mass and exact center of mass 
energy is not known. Therefore, it is sensible to concentrate on Lorentz invariant quantities like 
the invariant mass. The invariant mass of n particles with energy Ei and 3-momentum pi is 
frame independent, and given by 

m 2 = (x>) 2 -(x>y (5.1) 

The lowest limit of M is in any case the sum of all n particle rest frame masses. If there is no 
external force the invariant mass of a system of particles does not change, even if some particles 
decay or annihilate. In the following calculations the speed of light c and Planck's constant h 
are set to one. 

Invariant mass distributions carry information about the particles at the end of a decay chain 
as well as about the decaying particles. For this study, the kinematic formulae and configurations 
presented in this chapter are relevant. In general it is necessary to distinguish between massive 
and massless particles. If the energy of a particle is a few magnitudes higher than its mass, 
it is a good estimation to treat it as massless. Therefore, in the following calculations [20] all 
Standard Model quarks and leptons - except the top quark in section 15.31 - are chosen to be 
massless. In the Monte Carlo simulations, however, the correct masses are taken into account. 

5.1 Sequential two-body decays 

Only two-body decays have to be taken into account for the processes in this study. In the 
two-body decay of a particle M into mi and rri2, the energy, E* , and 3-momentum, p*, in the 
rest frame of M, which is denoted with a " * " , are given by 

M = E* 1 + E*,\p* 1 \ = \p*\=p*. (5.2) 

Squaring the first equation gives 

2E\El = M 2 -m\-m\- 2p* 2 . 
Squaring again and extracting p* leads to 

p* 2 = ^K M2 -™?-™i) 2 - 4 ™M 

[M 4 - 2M 2 (m\ + ml) + (m 2 - ml) 2 } 



AM 2 



[M 2 - (mi +m 2 ) 2 ][M 2 - (m ± - m 2 ) 2 ]. 



4M 2 
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The energy becomes 



77^*2 *2 i 2 

E l = p + mi 



[M 4 - 2M 2 (m? + ml) + (m? - ml) 2 + 4M 2 m 2 ] 



AM 2 



[M 2 + m\ - ml} 2 . 



AM 2 



In summary: 

M 2 + m| - m\ 

, ., , ., [(M 2 - (mi + m 2 ) 2 ) (M 2 - (mi - m 2 ) 2 )] 1/2 

IPil = IP2I = — 7T^ ( 5 - 4 ) 

The formulae are considerably simplified if one of the decay particles is massless. Assuming that 
m.2 = one obtains: 

_,„ M 2 + m 2 M 2 -m\ M 2 - m\ 

E\ = -^ , El = -^ , |p;| = |^| = -^. (5.5) 

The transformation from the rest frame of M to the lab system is defined by two kinematic 
variables of M, the velocity (3 = q/E and 7 = E/M, where q and E are the lab frame momentum 
and energy of M. The Lorentz transformation to the lab system is most easily decomposed into 
a longitudinal and a transverse part: 

Ei = j[E* + Pp*cos6*] 
Pn = Pri = P*sin9* 
p Li = 7 l(3E*+p*cos6*} (5.6) 

with #2 = 6\ + 71 "- The maximum and minimum values of Ei and pn are obtained for cosO* = +1 
and —1 respectively. 

The investigated leptonic X2 decay chain in the analysis presented here is 

X^^ + P ,l ± ^xi + l ± . (5.7) 

The calculations are valid for taus as well as for other leptons /, which are assumed to be massless. 
Moreover, maxima and minima of effective masses correspond to collinear configurations, where 
the particles are emitted along or opposite to the direction of the Lorentz boost. In this case, 

E n = lE* n (1 ± {3) . (5.8) 

As the Lorentz transformation parameters are 

M 2 - M? M 2 + M? 



P ' ™ 2 +M? ' 7 2i\L>Mr ' (5 ' 9) 

X2 l Xi 1 



this yields to 



-y(l±8) = (M 2 + M?) ± (M 2 - M?) (5.10) 

so that the maximum energies associated with the boost direction and minimum energies asso- 
ciated with the opposite-to-the-boost direction are: 

M-o . Mr 

E% ax = ~^E\ 2 , est = xr-^2- (5.11) 

Mr M^o 
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Similarly, for the momentum with the sign being measured in the boost direction it is 

Pl2 = 1 Et 2 (f3±l) (5.12) 

leading to 

Pl2 = h l2 > P\2 = ~ h l2 • (b.16) 

These expressions considerably simplify the later calculations. In the following the kinematic 
limit for both leptons is calculated. The effective mass can be computed in the rest frame of X.2- 
The z-axis is chosen along the direction of the slepton, as depicted on figure IS~T1 



~0 






~0 

12 

Figure 5.1: Kinematics of the \2 decay via a slepton I. 

The dileptonic effective mass is given by 

Ml = (E n + E l2 f - (- m + PLl2 f - (0 + P ti2? (5.14) 

= 2E lx E l2 + 2 Pll p Ll2 = 2En(Ei2 + Vhn) 

where the indices L and T denote the longitudinal and transverse components. As particle 11 is 
massless, the energy and momentum of 11 and I are given by the expressions (|5.5|) and similarly 
for particle 12 in the rest frame of I. 

Mf - M\ 
Eh = \ph\ = -^^ (5-15) 

The Lorentz transformation of the massless particle 12 into the rest frame of X2 becomes 

E l2 = yEf 2 (l +j3cos6*) 
PL12 = -fEf 2 (P + cos9*) (5.16) 



from which 



The effective mass then takes the simple expression 



E n +PL12 = 7^*2(1 + /?)(! + cos9*) 

M 2 
= ^(M?-M 2 )(l + cos9*)- (5.17) 



M 2 - Mf M 2 

Ml = 2^2— — -^(M? -M 2 )(l + cos6* 

11 2M*o 2M r y l *i A 

r2 T\/r2\(*/r2 fl/f 2 



{M'-M){M-M') 

— ^ — =-* 1^(1 + cose*) (5.18) 



2Mf 
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with the upper endpoint of the mass distribution given by 



MT X = M n 



\ 




M xf 



M x\ 



(5.19) 



This corresponds to the configuration where the two leptons are emitted back-to-back in the 
rest frame of the X.2- Since the I is a particle with spin it should decay isotropically. Thus, 
formula (J5.18J) shows that the distribution in Mu increases linearly with Mu leading to a sharp 
edge at the kinematic limit given by formula (|5.19|) . 

The following section shows a summary of all formulae for the endpoints which are used in 
the analysis and the related configurations. They are calculated in the same way as described 
above. 



5.2 Formulae and kinematic configurations 

Here the endpoints available for the decay chain q — > q + X.2 , X2 — > 11 + I , I — ► 12 + Xi are 
summarized. 



1. MS"" 



M max 



M X% 



\ 




(5.20) 



h 



X 2 



I 



Xi 



I 



2. M^ x (1st lepton) 



M hq 



M- q 



A 




M? 



Mf 



(5.21) 



q 



1 

I 



X 2 






h 
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3. MJ^j* (2nd lepton) 



M^ = Mg 



\ 



M xf 

Mg 2 



M xf 



(5.22) 



q q Xi 



k 



Xi 



4. MJf q ax , for M 2 < MgM^o 



MJ%*=M q 



\ 



M x% 

Mg 2 



M xl 

M xl, 



(5.23) 



q q 



I 



Xi 



h 



Xi 



I. 



5. (M Uq + Mb,)" 



^2 



(5.24) 



q q 



Xi 



I 



k 



Xi 



I. 
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Four of these five endpoints are in principle necessary to determine all involved masses, provided 
that several consecutive decay channels are open (long decay chains). However, it will not always 
be possible to measure them precisely. Therefore, it will be sensible to use all available endpoints. 

5.3 Endpoints for top quark events 

The mass of the top quark is very large and therefore not negligible. Thus, not only upper 
endpoints but also lower endpoints arise in the invariant mass distributions. In the analysis 
presented in section 16.2.131 to 16.2.151 the top quark is involved. The corresponding decay chain 

is 

i^t + X2^>t + l + l->t + l + l + Xi. (5.25) 



With 



M 2 + m? - Mi 



-j. i ■ ■ t'f 1VJ - ~o i 

E t - — —r , Pt = \M 2 " m 2 t (5.26) 



2M- t 



M? + M 2 - m 2 t 

E & = ^T , PvS = Mo ~ Wo (5-27) 



the endpoint formulae are the following: 



M |o m? M? E^o + p^o 

V m ' * tV M? M? IK M 2 J Mr y ' 

t t X 2 t 



x\[ I + ,''1 



<i >2 ) 



1 (■, ^ , m i\2 



(M$ n ) 2 = m 2 t + Mhl - -%- + ^){l - — L-)-*L—*L (5.29) 

K lit J t t y M 2 M 2J\ M 2 ) M V > 

t t xl % 



M |o m 2 M l° #v°+Pv° 

{MJSr? = m 2 + M?(l-^ + -±)(l--^) *\ g * 2 (5.30) 



Mi Mf JK M 2 ' Mr 
tt i <■ 



xi| | + 



W) 



M l° m 2 M? M 2 E+o-p+o 

(M,™ m ) 2 = m 2 t +M?(l %. + !!*) * h ^±)^2 — Zx* (5.31) 

v l2t ' * * v M? M? J M 2 K M 2 ' Mr y ' 

t t x 2 1 % 



"5" 



W) 



M'i m 2 M 2 E-o + v-q 



(M^ ax ) 2 = m f 2 + M^(l-— ^ + — 1)(1-— ^) *' X2 (5.32) 

v "* ; * * v M? M?' v M 2 J Mr y ' 

t t xl % 



M 2 m 2 M 2 E-o - p-o 

(M^) 2 = ml + M|(l - -^ + -|)(1 - -|i) *L_ ^ (5.33) 
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The configurations for the upper endpoints here are the same as for the massless case described 
in section HT21 The lower endpoints are realised with a configuration where all detected particles 
go into the same direction. 
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Chapter 6 

Searches for the X2 



Since mSUGRA is a Supersymmetry breaking scenario within the minimal supersymmetric 
extension of the Standard Model favoured by many theorists, it is of high importance to develop 
analysis methods for this model. Although the measurements at LEP (see section H~T]) and with 
the WMAP satellite (see section 14. 2[) give a strong constraint on its parameter space there is 
still a huge number of open possibilities with very different phenomenologies. Many studies for 
mSUGRA have been performed in the past years. Most of them concentrate on methods for 
models with tan/3 lower than 10. Actually, there is no reason for preferring low tan/3-models. 
It has been pointed out in several recent publications that due to the different phenomenology 
at different tan/3 values, studies for high tan/3 values are very important. 

Since there is no evidence for a heavy charged stable particle in the universe and assuming Im- 
parity conservation, the lightest supersymmetric particle in mSUGRA is the stable neutralino X\. 
A neutral weakly interacting stable particle escapes the detector, thus, a direct mass reconstruc- 
tion is not possible. However, it is useful to search for the next-to-lightest neutralino, the X2, 
which decays in detectable particles within the detector. It is mainly produced in gluino and 
squark decays. The focus here is on the leptonic cascade decay channel shown in figure IfTTl 




Figure 6.1: Leptonic cascade decay of %z. 



It dominates in areas of the parameter space where M^o is larger than Mj. Results from 

studies for X2-searches in mSUGRA at low tan/3 [21] are presented very briefly in the following. 
Then the results of the Monte Carlo analysis at tan/3 = 35 performed in this work are discussed. 
The latter shows that in principle there are ways to reconstruct the involved sparticle masses in 
that area of the parameter space. 

6.1 Searches at low tan/3 

Leptonic decays of the X2 have a useful kinematic characteristic. The dilepton invariant mass 
spectrum has a sharp edge near the kinematic upper limit which is given by formula 1)5.20(1 . 
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Thus, constraints on the involved sparticle masses which are independent of the underlying 
model can be obtained by measuring the endpoint of this distribution. The most straightforward 
signature for selecting the X2 decays is provided by the topology with two same-flavour opposite- 
sign leptons accompanied by large missing transverse energy which comes from the Xi. It is the 
advantage of low tan [3 models that many of the charged leptons coming from the X2 are electrons 
and muons which can be fully reconstructed. Figure E21 shows two representative examples of 
X2 branching ratios as a function of m 1 / 2 within the WMAP constrained parameter space. 
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Figure 6.2: Dominant branching ratios along the WMAP line for the X2 as a function of mi/2 
at tan/3 = 5 (left) and tan/3 = 10 (right), /i is positive and Aq is set to zero. The locations 



of the updated benchmark points A' (mo = 107, m 



1/2 



600), B'(57,250) and C'(80,400) are 



indicated [5]. Here leptons 1 are electrons and muons whereas taus are shown separately. 



In the two opposite-sign leptons + EJp tss + (jets) channel the largest Standard Model 
background is ti production, with t — > bW, and both Ws decaying into leptons, or one of the 
leptons coming from a W decay and the other from the 6-decay of the same t-quark. Defining 
the significance a as a = (Nev ~ Nb)/\/~Nev it has been shown that the observability of the 
kinematic endpoint varies from a = 77 at mo = 200 and m\/2 = 160 to a = 27 at mo = 60 
and mi/2 = 230 for an integrated luminosity Li nt = 10 3 pb~ l [21]. For this luminosity, the 
appearance of the edges in the distributions is already sufficiently pronounced in a significant 
part of the mo- 771 i/2"Pl ane - The position of the kinematic edge can be measured with a precision 
of about 0.5 GeV. Figure lfi~31 shows different realisations of such an edge very clearly. 
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Figure 6.3: The solid line shows the invariant mass distribution of two same-flavour, opposite- 
sign leptons of X2 decays at various (mo^i/2) points for tan/3 = 2 and \i < with 
Li n t = 10 3 pb^ 1 . The Standard Model background is given by the dashed line. 



With increasing mo and mi/2 the cross sections are decreasing and therefore a higher in- 
tegrated luminosity is needed. In order to perform a full mass reconstruction of the sparticles 
further analysis methods are necessary which are described to some extent in [21]. 



6.2 Searches at tan/3 = 35 



of) 



The kinematic characteristic of the leptonic X2 decay does not change with increasing tan/?. 
However, the branching ratios do. Simulations made with ISASUGRA show that the decay 
into taus became more and more important with increasing tan/3. Figure l6~H illustrates the X2 
branching ratios, evaluated at tan/3 = 35. 
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Figure 6.4: Dominant branching ratios along the WMAP line for X2 as a function of m.1/2 at 
tan/3 = 35, \i > and Aq = 0. The location of updated benchmark points I'(mo = 181,mi/2 = 
350) and J'(299,750) are indicated [5]. Here leptons 1 are electrons and muons whereas taus are 
shown separately. 



Neither with CMS nor with ATLAS will it be possible to fully reconstruct taus since they 
decay into at least one u T which escapes the detector. Therefore, the kinematic endpoint in the 
dilepton invariant mass distribution will in general be smeared out. It is the aim of the following 
study to investigate the X2 discovery potential in this region of the mSUGRA parameter space. 
Furthermore, the goal is to use the X2 discovery for sparticle mass reconstructions. 



6.2.1 Choice of parameters 



At the moment it is not known whether the MSSM with mSUGRA describing the Supersym- 
metry breaking is the right extension of the Standard Model. However, there are experimental 
(see section l4~Tl and EO|) and theoretical constraints on the mSUGRA parameter space. For the 
development of analysis methods it is necessary to choose parameters in a characteristic region 
which is not yet excluded by any experiment and which is accessible for future experiments like 
the Large Hadron Collider. 
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Figure 6.5: Points for the analysis in the mQ-rrii / 2 plane in an area around the benchmark point 
F. Fhe upper left side is excluded due to cosmological constraints. 



Fhe parameters in this analysis are chosen with the following motivation: 

1. tan/3 = 35: Fhere is no experimental limit which excludes high tan/3 models. It is still 
necessary to develop methods for the X2 discovery at LHC assuming these models. It is 
especially necessary to determine whether it is possible to discover the X2 in the r-channel 
which is favoured at tan/3 = 35. 

2. rriQ £ [167, 193], m\ / 2 £ [340,360] around F in figure l6~H In that region the branching 
ratio BR(X2 — ► f{ r^) is with over 95% especially high. For the analysis 25 points in that 
region are chosen. In order to consider the WMAP results and the cosmologically excluded 
region the points are not equidistant which can be seen in figure l6~5l 

3. Aq = 0: A non zero value would only be necessary for the electroweak symmetry breaking 
if the top mass was lower than the measured value. However, it would lead to different 
sparticle masses. 

4. /i > 0: Fhere is no special reason to choose /i > 0, but the results are not very sensitive 
to that choice [20]. 

Fhe masses of the sparticles relevant for this study in the chosen area of the mSUGRA parameter 
space are given in table HTT1 
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Masses in GeV for (mo,mi/2) 


Particle 


(167,340) 


(187,340) (181,350) 


(173,360) 


(193,360) 


9 


812.06 


812.49 834.16 


855.61 


856.23 


d\/s\ 


759.81 


764.82 782.61 


800.09 


804.65 


(I2/S2 


730.15 


735.04 752.02 


768.81 


773.64 


ui/ci 


755.24 


760.29 778.17 


795.75 


800.34 


U2/C2 


731.79 


736.67 753.80 


770.73 


775.06 


h 


657.94 


661.48 678.14 


694.20 


697.50 


b 2 


704.72 


708.19 724.93 


740.96 


744.21 


h 


565.70 


568.46 583.27 


598.06 


600.47 


t 2 


729.25 


731.74 747.90 


763.69 


766.14 


x\ 


257.10 


257.23 265.52 


273.74 


273.92 


Xx 


133.93 


133.99 138.11 


142.28 


142.41 


h 


135.74 


154.52 150.14 


144.24 


162.56 



Table 6.1: Sparticle masses in GeV calculated with ISASUGRA 7.69. Aq is set to zero, \x is 
positive and tan/3 is 35. 



6.2.2 Production, decay and selection of X2 

There are different possibilities of producing X2, and there are different channels in which the X2 
can decay. The four X2 decays with the highest branching ratios are shown in table l6~2l Tabic lo~3*l 
shows eight processes with the highest cross sections which either produce a neutralino directly, 
or indirectly in decays of squarks and gluinos. 





Figure 6.6: Cascade decay of g and q. 



The latter processes are used in this study since there are one or two additional quarks 
in the decay chain which carry information about the involved sparticle masses. They can be 
combined with both leptons in order to get more kinematic endpoints. The processes are shown 
in the two diagrams in figure ESI Events with directly produced X2 can only be used for finding 
the dilepton kinematic endpoint. The cross sections are calculated using PYTHIA 6.220 in 
combination with ISASUGRA 7.69 at the given (mo, mi/ 2) values. The reliability of these cross 
sections is described in appendix IA. 41 For the indirect production, the squark branching ratios 
of the decay into the X2 - shown in table l6~H - are important. The left handed squarks show much 
larger branching ratios for the decay into wino-like X2 ■ The gluinos will in most of the cases first 
decay into a squark and a quark and then through the squarks produce X.2- Table HT5l shows 
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that a significant amount of left squarks is produced through gluino production. Additionally, 
it can be seen that the sbottom channel is favoured in the gluino decay. 

In this study the analysed data samples correspond to 100 000 generated events in which the 
tau decay is specified. Depending on the values for mo and mi/ 2 this leads to between 24278 
and 25379 selected events. The event selection in this study is at the Monte Carlo level. In 
the real experiment it will be an additional challenge to select and reconstruct the taus and the 
associated quarks. Here an event is defined as follows: 



Two opposite-sign charged rho mesons (see section Rj. 2.4(1 which come from two opposite- 
sign tau leptons. These tau leptons must come from the cascade decay of the X2 which is 
produced through squarks. At this level the exact origin of each particle is known. In order 
to take the experimental reconstruction of the rho mesons into account, the distribution 
of the rho mass is included into the simulation. In the real experiment two opposite-sign 
isolated tau-like jets have to be searched with a 7r° and a ir^ 1 in the detector having an 
invariant mass around the p mass. 



2. The quarks used in the analysis come only from squarks which decay into a X.2- The 
flavour is known and also the precise 4-momentum of the quark is used for the analysis. It 
is mostly called "associated" quark in this study. In the experiment at least one additional 
jet is needed. 



3. No cut on the E™ lss 



has been applied. However, in the real experiment E™ lss should 
required to be greater than 100 GeV. This is the value chosen in most of the low tan/3 
studies. It has to be investigated whether this limit is also a good value for high tan (3. 



Branching Ratio in % for (too^i^) 


Decay 


(167,340) (187,340) (181,350) 


(173,360) 


(193,360) 


X2 -► ff + TT 


96.8 96.2 96.0 


96.0 


95.2 


X 2 '->X U i + h° 


2.1 2.6 2.8 


2.9 


3.6 


x^xU z° 


0.7 0.9 0.8 


0.7 


0.9 


X2 -► 1% + P 


0.4 0.2 0.3 


0.4 


0.3 



Table 6.2: The largest branching ratios of the X2 for the chosen mSUGRA parameters. Aq is 
set to zero, [i is positive and tan/3 is 35. The letter I denotes electrons and muons. 
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Cross Section 


(pb) for (m 


,m 1/2 ) 






Subprocess 


(167,340) 


(187,340) 


(181,350) 


(173,360) 


(193,360) 


9+9 


5 + 5 


0.636 


0.634 


0.526 


0.434 


0.435 


Qj+9 -► 


<?jL + 3 


1.865 


1.824 


1.578 


1.366 


1.338 


Qj + 9 -» 


9ji? + 9 


2.030 


1.987 


1.718 


1.490 


1.470 


g + <? -»• 


xl + xt 


0.271 


0.276 


0.241 


0.214 


0.217 


% + <?j -»• 


<7iL + 9jL 


0.502 


0.496 


0.443 


0.399 


0.393 


Qi + 9j -► 


QiR + <7j.R 


0.585 


0.575 


0.515 


0.466 


0.459 


% + Qj -► 


(&L + Qjr) 


0.467 


0.454 


0.398 


0.365 


0.356 


+ 


{QiR + QjL) 












% + <7j -► 


{qiL + qjR.) 


0.518 


0.509 


0.444 


0.398 


0.384 


+ 


(qiR + qjl) 













>-'.! 



Table 6.3: Most important processes for direct and indirect (in cascade decays of g and q) X2 
production at LHC. Aq is set to zero, /j, is positive and tan/? is 35. 





Branching Ratio in 


% for (rao,rai/ 2 ) 






Decay 


(167,340) (187,340) 


(181,350) 


(173,360) 


(193,360) 


d L ^X2 + d 


31.1 31.0 


31.1 


31.2 


31.2 


d R ^X°2 + d 


0.3 0.3 


0.3 


0.2 


0.2 


UL —> X2 + U 


31.5 31.5 


31.6 


31.6 


31.6 


UR — ► X2 + U 


0.3 0.3 


0.3 


0.2 


0.2 


SL —> X2 + S 


31.1 31.0 


31.1 


31.2 


31.2 


SR — >■ X2 + S 


0.3 0.3 


0.3 


0.2 


0.2 


CL —> X2 + C 


31.5 31.5 


31.6 


31.6 


31.6 


CR -*■ X2 + C 


0.3 0.3 


0.3 


0.2 


0.2 


&i -► X2 + 6 


29.3 29.0 


28.5 


28.1 


27.9 


&2 -► X°2 + b 


7.3 7.0 


7.0 


7.1 


6.9 


fl -► X2 1 + 1 


14.3 14.3 


14.3 


14.4 


14.3 


f 2 -» *2 + t 


9.1 9.2 


9.3 


9.4 


9.5 



•r,o 



Table 6.4: Squark branching ratios for the decay into X.2- Aq is set to zero, // is positive and 
tan/3 is 35. 
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Branching Ratio in 


% for (m ,mi/ 2 ) 






Decay 


(167,340) (187,340) 


(181,350) 


(173,360) 


(193,360) 


g-^X 2 + t + t 


0.01 0.01 


0.01 


0.01 


0.01 


9 -> d L + d 


3.5 3.2 


3.3 


3.5 


3.2 


g -»• d R + d 


8.1 7.9 


8.0 


8.0 


7.8 


g -^u L + u 


4.1 3.8 


3.9 


4.0 


3.8 


9 ->• u R + u 


7.8 7.6 


7.7 


7.7 


7.5 


9 ~* sl + s 


3.5 3.2 


3.3 


3.5 


3.2 


9 -»• s R + s 


8.1 7.9 


8.0 


8.0 


7.8 


9 -> c L + c 


4.1 3.8 


3.9 


4.0 


3.8 


9 -*• cr + c 


7.8 7.6 


7.7 


7.7 


7.5 


5 -» h + 6 


24.8 25.8 


24.8 


24.0 


24.8 


g ^b 2 + b 


13.9 14.2 


14.0 


13.8 


14.1 


g-^h + t 


14.3 15.0 


15.5 


15.8 


16.5 


~g^t 2 + t 


0.0 0.0 


0.0 


0.0 


0.0 



Table 6.5: Gluino branching ratios for the given {mo,mi/ 2 ) values. Aq is set to zero, // is positive 
and tan (3 is 35. For the chosen parameters the g cannot decay into i 2 +t because M g < Mj 2 +M t . 



6.2.3 Methods to find the endpoint 



c0 



Due to the presence of neutrinos in the investigated X 2 decay, the distribution of the invariant 
mass of all visible X 2 decay products is fundamentally different from the case where the X 2 
decays into a X\ and a pair of electrons or muons. The neutrinos smear out the kinematic 
endpoint, which is shown in figure I7T71 Therefore, it is necessary to define a method how to find 
the kinematic endpoints in the measured invariant mass (Mjjvv) distributions. In the following 
two "simple" types are presented. "Simple" means here that they do in principle not depend on 
the model. However, the included corrections for the improvement of the measurement and the 
obtained systematic uncertainties are based on mSUGRA Monte Carlo simulations. 

Linear fit 

First, there is the possibility to perform a linear fit f(MjNv) to the invariant mass distribution 
near the kinematic limit. The invariant mass range in which the fit is performed needs to be 
chosen carefully in order to obtain a sensible fit. Certainly, the region is between the bin with the 
highest content and the end of the tail. Good results have been achieved using an area starting 
at the bin with 50% and ending at the bin with 5% of the maximum bin content. Definitely, 
this definition is only meaningful if no background is included. The measured endpoint Epf 1 
for each point i of the ./V = 25 investigated points in the parameter space is then defined as the 
intercept point of the fit line and the abscissa. A{ is a free parameter which is not needed for 
the endpoint measurement. 



f(M INV ) = A i -(M INV -B i ) 



Ep: 



Al 



B; 



(6.1) 

(6.2) 



The mean statistical uncertainty on Ep is 



o~statEp 



N 



(6.3) 
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with 5Bi being the statistical uncertainty of Bi. The systematic shift S of the endpoint measure- 
ment is given by the mean difference of the theoretical endpoint Epf and Epf 1 , which includes 
a systematic and a statistical uncertainty 



££i Epf - Epf 

N 



(6.4) 



68 = f^^^-Sf (6B) 

The corrected measured endpoint is then 

Epf' c = Epf + S. (6.6) 

The pure systematic uncertainty on the endpoint measurement is then 



8 syst Epf = ^(5S) 2 -(5 stat Epf) 2 (6.7) 

Certainly, the calculated S can and will change when a different region in the mSUGRA param- 
eter space is chosen. 

Gaussian fit 

Second, the shape of the distribution can be investigated. In the chosen region the neutrinos 
smear out the visible invariant mass distribution near the kinematic limit. As a result the 
maximum of the invariant mass distribution is at a mass much lower than the endpoint. The 
measurable maximum Gi is here defined as the maximum of an iterated gaussian fit. The 
maximum position is determined as follows: 

1. The bin with the maximal bin content is the centre for the first gaussian fit. Unless 
otherwise noted for the investigated invariant mass distributions a sensible area for the 
gaussian fit is a symmetric area around the centre of 25% of the centre value. 

2. The resulting gaussian maximum is then the input for the centre of the second fit. 

3. It turned out that three iterations are enough for a stable maximum. The last (out of 
three) gaussian maximum is taken for Gi. 

Gi has a statistical uncertainty 5Gi. If there is a correlation between the maximum and the 
endpoint which is independent of itiq and mi^ over a large area, then it is possible to get the 
endpoint value out of the maximum measurement. In the following it is assumed that there 
exists a linear dependence between Gi and Epf. A linear least squares fit of Epf as a function 
of Gi ± SGi with the fit function f(G) = C ■ G + D gives the values for this dependence. The 
measured endpoint Epf with its mean statistical uncertainty 5 s t a tEp is then 

Epf = C-Gi + D, (6.8) 

6 stat Ep M = C- ^ 5Gi . (6.9) 

The corrected endpoint, the systematic shift with its uncertainty and the pure systematic un- 
certainty are equally defined as in the case of the linear fit. Since the gaussian maximum should 
be zero if Epf is zero it is interesting to check whether D is close to zero or not. This value 
gives a hint for the stability of the assumed correlation on further changes in mo and mi/2- I n 
the following analysis therefore also the slope C and the resulting uncertainties on the endpoint 
measurement are given if D is forced to zero. It is obvious that the gaussian method relies on 
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Monte Carlo simulations in order to find C and D whereas the measurements of a sharp edge 
at low tan/3 are model independent. 

The difference between Epf and Epf + S in the linear and in the gaussian case may have 
in some cases a mo and m^ft dependence. This is discussed for each invariant mass distribution 
separately. For the analysis the binning and the choice of the fit area have been optimised such 
that they are most universally valid within the investigated area of the mSUGRA parameter 
space. 

6.2.4 Selection of the leptonic decay 
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Figure 6.7: The red line shows the visible dilepton invariant mass distribution in the /r^-channel, 
the green line in the 7r -channel and the blue line in the \x -channel for mo = 181 and TO1/2 = 
350. The gray (1) and the black (2) line show the distribution of the 3-prong channel. The first 
includes the r ± — ► is^(Pv T and the direct r ± — ► tt^tt^tt^Vt- decay whereas in the second only 
the latter is included. 



Unlike electrons and muons, the tau leptons cannot be fully reconstructed, since they decay 
in at least one neutrino within the detector. However, it is the aim to find the invariant mass 
upper limit of both tau leptons in the leptonic X2 decay. The kinematic endpoint of the visible 
dilepton invariant mass distribution is the same as the limit of the true invariant mass because 
neutrinos are massless. However, its shape is strongly dependent on the chosen visible particles. 
The branching ratios for the tau decays in detectable particles are shown in table WM For the 
r — ► p v r channel, the visible mass of the X2 decay products tends to be higher than for the 
other channels, except the 3-prong, which can be seen in figure l6~71 There the neutrinos get less 
energy than in the other decay chains because of the relatively large rho mass. The branching 
ratio of 25.41%, compared to the 14.57% of the 3-prong channel, is a further advantage of this 
channel. Therefore, in the following analysis the r^ — ► p^v T channel is chosen. The /9 ± used 
here are at Monte Carlo level with a mean mass of 766 MeV and a width of about 126 MeV 
given by PYTHIA. 
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T 



± 



p ± u T 






e^v e v T 



t —> h h h + > neutrals v T 
(3-prong, ex. Kg — » 7r + 7r~) 



T — > 

Other 



7T i/ T 



(25.41%) 
(17.36%) 
(17.84%) 
(14.57%) 



(11.06%) 
(13.76%) 



Table 6.6: The branching ratios of the tau lepton [12]. 



6.2.5 The kinematic limit of M p ± pT 
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~0 
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~0 




T 2 



The mass of the f is quite similar to that of the Xi for small values of mo in the chosen 
mSUGRA parameter space. Formula (J5.20J) shows that in this case the endpoint is close to zero. 
Therefore, the theoretical values for the dilepton kinematic limit given in figure IB~%1 have a strong 
mo and m\/2 dependence. In the decay chain of the X2 at high tan/3 there will in about 6% 
of the cases be two opposite-sign rho mesons which can be used for the endpoint measurement. 
The shape of the invariant mass spectrum of both rho mesons shows that there is no sharp edge 
near the kinematic limit since it is smeared out due to the neutrinos. However, there are two 
possibilities to measure the endpoints (see section r6.2.3[) . In the following analysis the whole 
data samples with about 25000 events are used to determine the corrections and systematic 
uncertainties for both methods. 

First, it is possible to perform a linear fit to the tail of the distribution. An example is 
given in figure E3 The measured endpoint mass is too small for every point in the investigated 
region. The mean statistical uncertainty 5 s t a tEp on this measurement for these large data 
samples is 0.2 GeV and therefore negligible. The mean correction S is 9.0 GeV which can be 
seen in figure IfTTUl Indeed, the uncertainty 5S = 1.8 GeV on this correction is small, but as can 
be seen in figure l6~TT1 the real needed correction should increase with increasing tuq. The pure 
mean systematic uncertainty 8 sys tEp M is 1.7 GeV. 

The second method can be applied for this distribution since the gaussian fit is stable for all 
25 points and the assumed linear dependence between G{ and EpJ is sensible (see figure l6~T2Jl . 
If the linear fit of EpJ as a function of G\ ± 5Gi is forced to go through the origin, the resulting 
slope C is 3.1. The statistical uncertainty S s t a tEp on the endpoint measurement has then a 
mean value of 2.2 GeV. However, in this case a systematic shift S of 1.5 ± 2.3 GeV is needed 
in the whole region to improve the endpoint measurement (see figure I6.13|) . If the fit is not 
forced to intersect the origin the value C changes to 3.4, no systematic shift S is needed and 
5S decreases to 1.9 GeV. The fit here intersects the ordinate at D = —6.1 GeV. The mean 
statistical uncertainty 5 s t a tEp on the endpoint measurement increases slightly to 2.4 GeV. 
The mean pure systematic uncertainty is not defined here because SS is dominated by the 
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statistical uncertainty in the data samples. The difference of the theoretical endpoint EpJ and 
Epf 1 shown in figure KTTH has no significant too and mi/2 dependence. 



Value 



6 s tatEp 

S 

5S 



w 



Linear Gaussian FO Gaussian NFO 



0.2 
9.1 
1.8 



2.2 
1.5 
2.3 



2.4 



1.9 



Table 6.7: A summary of the calculated mean uncertainties and systematic shifts in GeV. 
5 s tatEp is the statistical uncertainty on the endpoint measurement, S is the mean difference 



of the measured and the theoretical endpoint with the root-mean-square deviation 5S defined 
in formula 16.41 and 16.51 FO denotes the results with fixed origin and NFO with non- fixed origin 
in the Gi-Epf plane. 
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Figure 6.8: Theoretical kinematic endpoints in GeV for the invariant mass of both opposite-sign 
rhos coming from the X2 and the f respectively. The endpoint masses decrease significantly for 
decreasing values of toq . 
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Figure 6.9: Example for the linear and the gaussian fit at m® = 181 and mi/2 = 350. In this 
distribution 25071 events are represented. The triangle shows the theoretical kinematic endpoint 
at 85.8 GeV. 
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Figure 6.10: The difference of the theoretical and the measured M p ± p ^ endpoint in GeV obtained 
with the linear fit at the 25 investigated points in the mo-mi/2 plane. It shows a mean value of 
9.1 GeV with a root-mean-square deviation of 1.8 GeV. 
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Figure 6.11: The theoretical M p ± p ^ endpoint minus the measured endpoint in GeV obtained with 



the linear fit after applying the constant correction of 9.1 GeV. The values decrease significantly 
for decreasing values of mo- 
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Figure 6.12: The theoretical dilepton kinematic endpoints as a function of the measured gaussian 
maximum for all 25 investigated points in the mo-mi/2 plane. The gray line is forced to go 
through the origin whereas the black has an optimised ordinate value. 
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Figure 6.13: The difference of the theoretical and the measured M p ± p ^ endpoint in GeV obtained 
with the gaussian fixed origin method. The measured values are on average lower by 1.5 GeV. 
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Figure 6.14: The difference of the theoretical and the measured M p ± p ^ endpoint in GeV obtained 
with the gaussian non- fixed origin method. 
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6.2.6 Selection of quarks 

In order to get more equations for the reconstruction of the sparticle masses it is necessary 
to include the quark which comes from the process qi — ► X2 + q%- Since the squark masses are 
flavour dependent the kinematic limits are as well. In the detector it may in principle be possible 
to distinguish up, down, charmed and strange quarks from bottom and top quarks. The top 
quark decays mainly into a bottom quark and a W boson. Therefore, by using the b-tagging 
it is in principle possible to divide the processes into two groups, light and heavy quarks. The 
upper kinematic limits of the invariant mass distributions with quarks depend on the most 
heavy squark. In the investigated region the heaviest squarks among the u, d, c and s are the 
left handed d\ and s\. Therefore only this mass can be reconstructed. The kinematic limits for 
events with cfe, S2, Ui,2 and ci 7 2 are hidden within the measured distributions. However, due 
to the proton structure the u\ dominates the distributions as can be seen in figure IHTTKl The 
heavy quark distributions are strongly dominated by the b\ (t\), however, 62 (£2) are heavier 
and determine the value for the upper kinematic limit. This problem is discussed in more 
detail in section 16.2.101 In this study bottom and top quarks are treated separately in order to 
understand their different contributions to the invariant mass distributions. This knowledge can 
then be used for more realistic simulations. The following results, which include information 
from quarks, are strongly dependent on the precision of the quark energy reconstruction. For 
the analysis quarks are used at the parton level. Indeed, this is a very crude approach but, 
nevertheless, the methods developed here may also be valid for the real experiment. 



6.2.7 The kinematic limit of M p ± pT „ for light quarks 
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The invariant mass distribution of the quark which comes from the decay q — ► q + X2 and the 
same two rho mesons used in section Rj . 2 . 51 carries further information which can be used for the 
reconstruction of sparticle masses. Figure E2D shows that the theoretical values of the kinematic 
limits have a very small ttiq dependence and a slightly stronger mj/j dependence in that area 
compared to the dilepton case. The reason is that the last factor in formula l|5.2HJI in that region 
is never close to zero in contrast to that in formula (|5.2(J[) . The uiq dependence is only caused 
by the change in the squark mass whereas the mi/2 dependence takes also the changes of both 
neutralino masses into account. For the following analysis data samples of about 18000 events 
have been used. They include the same events as used in the previous section with the constraint 
that the associated quarks have the flavour up, down, charmed or strange. An estimation of 
the kinematic limit can be achieved with a linear fit. An example is given in figure l6~T7l The 
measured value is too low for every point in the parameter space. As can be seen in figure IfTTKl 
and figure I6~T§1 the actually needed correction is not constant within the investigated area. The 
gaussian method cannot be applied to this case, since there is no linear dependence between the 
maximum and the endpoint - which can be seen in figure l6~2"Ul 
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Figure 6.15: The invariant mass distributions for two opposite-sign rhos with the associated light 
quark q (q — > q + X.2)- The gray triangle shows the theoretical endpoint for the u\ and c\ distri- 
bution at 624 GeV and the black triangle that for the d\ and s\ distribution at 628 GeV. The 
u\ distribution clearly dominates, but due to the small difference of both theoretical endpoints 
the effect is small, too. 



Value 


Linear 


SstatE P M 
S 

6S 


1.7 

14.0 

7.8 



Table 6.8: A summary of the calculated correction and the mean uncertainties in GeV for the 
endpoint measurement in the distribution of M p ± p ^ q . 
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Figure 6.16: Theoretical kinematic endpoints in GeV for the invariant mass of both opposite-sign 
rhos, coming from the X2 and the f respectively, and the associated quark from §x,di — ► s, d+X.2 
in the investigated area. 
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Figure 6.17: Example for the linear and the gaussian fit at m® = 181 and 77i 1 / 2 = 350. The 
invariant mass distribution is based on 18714 light quark events which are defined in section lfi.2.2l 
The triangle shows the theoretical kinematic endpoint of M p ± p ^ q at 627.8 GeV. 
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Figure 6.18: The theoretical M p ± p ^ q endpoint minus the measured endpoint in GeV for all 25 
points. The mean value is 14.0 GeV with a root-mean-square deviation of 7.8 GeV. 
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Figure 6.19: The theoretical endpoint minus the measured endpoint in GeV after applying a 
constant correction of 14.0 GeV. 
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Figure 6.20: The theoretical kinematic dilepton endpoints as a function of the measured gaussian 
maximum for all 25 investigated points in the rriQ-m\/2 plane. A unique mapping of the measured 
maximum to the theoretical endpoint can not be established since several theoretical endpoint 
values are associated with the same measured maximum. Additionally, some distributions are 
flat around the maximum which leads to large uncertainties. 



6.2.8 Two kinematic limits in M p ±„ for light quarks 
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As described in section I5.ll there are two different leptons in the investigated decay chain. 
The first comes directly from the X2 and the second from the slepton. Both leptons have 
different kinematic properties so that the endpoint of the invariant mass of first lepton with the 
light quark is not the same as that with the second lepton. However, in the experiment it will 
not be possible to know for each event which is the first lepton and which the second. In the 
following it is described how, nevertheless, both endpoints can be separated if the values are not 
too close to each other. For the latter constraint it is necessary to compare formula ()5.21|) with 
formula (|5.22|) . In the chosen area of m,Q and m 1 / 2 the f mass is much closer to the Xi mass 



45 



~0 



than to the X2 mass. Thus, the second endpoint is much lower than the first one and it will be 
possible to distinguish both endpoints. First, the method for finding the second endpoint will 
be described as its outcome is useful for measuring the first endpoint. 

For the measurements of the smaller second endpoint (given in figure I6.22J) the first rho 



gives a background which has to be avoided. In figure IH72T1 the invariant mass distribution of 
the single rho with one quark versus the invariant mass of both opposite-sign rhos is shown. 
On the left upper quadrant the events with the second rho dominate. Therefore, by ignoring 
the events with an invariant mass of both rho mesons being less than or equal to 55% of its 
upper limit the background can be well reduced. This cut value has been uniformly optimised 
for all 25 points. The ratio of the signal to the first lepton background between zero and the 
theoretical value for the second kinematic endpoint can thus be increased by a factor which is 
between 1.4 and 1.9, depending on the values for mo and m\/2- This improvement can be seen 
by comparing figure HT2"51 with figure l6~2"H For each individual point in the parameter space the 
signal to background ratio could even be improved by changing the cut value. For the following 
analysis data samples of between 5096 and 6806 events after the cut have been used. 
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Figure 6.21: The invariant mass of both opposite-sign rho mesons versus the invariant mass of a 
single rho with the associated quark in GeV. It is an example for the point I'. The right ellipse 
denotes the events with the first rho whereas the left ellipse includes the second rho events. The 
lower M p ± q endpoint is determined from the events with M p ± p ^ > 46 GeV (unshaded region). 



The linear fit method has to be adapted to this case. In order to find the endpoint on the 
background two linear fits are necessary. The first fit with the fit function /(-M/jvv) is defined 
the same way as in section lb. 2 . 31 except that the chosen region is from 95% to 45% with respect 
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to the maximum bin content. The second fit with the fit function g(Mj^v) iteratively searches 
the area right of the second lepton peak where the slope is bigger than —0.2 for a fit over seven 
bins. The minimal slope and the number of bins have been optimised for a stable second fit. 
The aim of this method is to find the second endpoint on the invariant mass distribution of the 



first rho with the associated quark. Here, the endpoint Epf is defined as the intercept point of 
both linear fits which are shown in figure 16.241 If the two linear fit functions are defined by 



f(M INV ) 

g(M INV ) 



A ■ M INV + B 
C ■ M INV + D 



the measured endpoint Epf is given by the intersection of / and g: 



Ep. 



The associated statistical uncertainty is: 



M 



Di-Bj 

Ai - Bt 



M\2 



(5 stat Ep M ) 



Sbo + ShB , (A-^) 2 -(^ A + ^ c ) 



+ 



{Ai — Cj) [A^ — L-i) 

(A-£i) 2 -<W*AB 



(6.10) 
(6.11) 



(6.12) 



(6.13) 



+2 



(A t - df 



The values 5 l XY are the entries of the error matrix which contains the correlated statistical 
uncertainties on A, B, C and D. Since g and / are not correlated the elements S l AC , 8\ D , S BC 
and S BD are zero. The result of this method is that the measured endpoint is for every point in 
the parameter space too low ( Figure lb. 25(1 . A correction S of 26.8 ± 8.0 GeV has to be added. 
The actually needed correction here has the same too dependence as that in M„± pT which can 
be seen in figure E2H1 

The gaussian method here remains the same except that the symmetric area for the fit 
changes from 25 % to 45 % of the centre bin value in order to achieve better results. The reason 
for this change is that the peak here has a smaller width. As can be seen in figure 16.271 it is 
reasonable to assume a linear dependence of Epf on G{ ± 6G{. If the origin is chosen to be fixed 
the slope C of the fit is 1.7. If the ordinate value is a free parameter the slope C decreases to 
1.5 with D = 22.7. The ratio between both slopes is 1.1. 

The calculated uncertainties and the systematic shifts for the linear and both gaussian meth- 
ods are given in table 16.91 . 



Value 



d s tatEp 

S 

5S 



w 



Linear Gaussian FO Gaussian NFO 



4.6 

26.8 

8.0 



3.0 

-3.9 

9.9 



2.7 
0.1 
2.9 



Table 6.9: A summary of the calculated mean uncertainties and systematic shifts in GeV for the 
second M p ± q 
section 16.2.31 



second M p ± q endpoint. FO is with fixed origin and NFO with non-fixed origin as described in 



When the second endpoint is determined it can be used to make a cut on the measured 
invariant mass in order to select solely first rho events. The uncertainty of that cut value has 
no effect on the following results since the tail of the distribution is not dependent on the cut. 
The only important aspect is that the fit algorithm for the linear fit takes the maximum of the 
first rho distribution and not the maximum of the second rho distribution. For the following 
analysis data samples of between 5142 and 16963 events after the cut have been used. The 
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large difference of the events in the data samples is caused by the very different ratio between 
both endpoint masses. However, the relevant number in the fitted region is in all data samples 
of the same magnitude. The theoretical first rho endpoints are shown in figure 16.301 In the 
investigated area the maximum of the first rho plus associated quark distribution is very close 
to the endpoint concerning the second rho. Therefore, it is only possible to make a linear fit 
at the tail of the distribution. An example is given in figure l6~3"Tl The systematic shift here is 
compared to the previous cases with 1.8 ± 2.4 GeV (Figure f6.32|) relatively small. Furthermore 
there is no significant mo and vn\ii dependence on that shift which can be seen in figure l6~3Hl 



Value 


Linear 


5statE P M 

s 

8S 


1.7 
1.8 
2.4 



Table 6.10: A summary of the calculated correction and the mean uncertainties in GeV for the 
first endpoint in M p ± q . 
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Figure 6.22: Theoretical kinematic endpoints in GeV for the invariant mass of the second rhos 
coming from the f^ and the associated light quark coming from the squark in the s\,d\ — ► 
s, d + X2 decay. The endpoints decrease significantly for decreasing values of ttiq. 
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Figure 6.23: Example for the invariant mass M ± Q distribution at itiq = 181 and m 
before the cut on M, 
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In this distribution 37428 events are represented. The triangle shows 



the theoretical kinematic endpoint of M ± at 288.0 GeV. 
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Figure 6.24: Example for the two linear fits / and g and the gaussian fit at ttiq = 181 and 
m i/2 = 350. In this distribution 6904 events after the cut (55%) on M„± „ T are represented. The 
triangle shows the theoretical kinematic endpoint of M ± at 288.0 GeV. 
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Figure 6.25: The difference of the theoretical and the measured M ± endpoint obtained with 
the linear fit in GeV for all 25 points. It shows a shift of 26.8 GeV with a root-mean-square 
deviation of 8.0 GeV. 
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Figure 6.26: The theoretical M ± endpoint minus the measured endpoint obtained with the 
linear fit in GeV after applying the constant correction of 26.8 GeV. The values decrease signif- 
icantly for decreasing values of ijiq. 
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Figure 6.27: The theoretical kinematic endpoints of M ± as a function of the measured gaussian 
maximum for all 25 investigated points in the mQ-mi / 2 plane. The gray line is forced to go 
through the origin whereas the black has an optimised ordinate value of 22.7 GeV. 
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Figure 6.28: The difference of the theoretical and the measured M ± endpoint in GeV obtained 
with the gaussian fixed origin method. The measured values are on average too large by 9.1 GeV 
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Figure 6.29: The difference of the theoretical and the measured M ± endpoint in GeV obtained 
with the gaussian non- fixed origin method. 
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Figure 6.30: Theoretical kinematic endpoints in GeV for the invariant mass of first rho coming 
from the X2 and the associated light quark coming from the squark in the Si,di — ► s,d + X2 
decay. 
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Figure 6.31: Example for the linear fit at ttiq = 181 and m 



1/2 



350 after applying a cut on the 



invariant mass M p ± a at the theoretical second endpoint value. In this example 9855 events are 



represented. The triangle shows the theoretical kinematic endpoint of M ± at 606.0 GeV. 
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Figure 6.32: The difference of the theoretical and the measured M± endpoint in GeV obtained 
with a linear fit for all 25 points. It shows a shift of 1.8 GeV with a root-mean-square deviation 
of 2.4 GeV. 
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Figure 6.33: The theoretical M ± endpoint minus the measured endpoint in GeV obtained 
with the linear fit after applying the constant correction of 1.8 GeV. The uncertainties show no 
significant tuq or mi/2 dependence. Thus, a constant correction is sensible in this case. 



6.2.9 The kinematic limit of M ± + M T for light quarks 
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If all endpoints could be measured very precisely four measured kinematic limits would be 
enough to reconstruct the four sparticle masses M^o, M^o, Mf and My-. But mainly due to 
the neutrinos even at parton level there are large statistical and systematic uncertainties for the 

H) depends on the previous four 



measured endpoints. Although this fifth endpoint (Figure | 
endpoints, its measurement improves significantly the uncertainties for the mass reconstruction 
which is performed in section 16.2.161 For the following analysis the same data sample as in 
section 15.2.71 is chosen. 

The linear fit (Figure d. 35j) measures the endpoints systematically too low by 51. 4± 12.8 GeV 
as can be seen in figure HU?6l The same mo dependence on the actually needed correction as has 
been determined for M p ± pT or M p ± pTq is shown in figure l6~3*71 

Figure IH.38I illustrates that the gaussian method is suitable to determine the endpoint. Even 
if the fit is forced to go through the origin, which leads to a slope value C = 1.9, good results 



54 



can be achieved. If the fit can optimise the ordinate value the slope C decreases to 1.3 with 
D = 254.3 GeV. The calculated values for the uncertainties and the systematic shifts are given 
in table IOTI 



Value 



Linear Gaussian FO Gaussian NFO 



w 



6 s tatEp 

S 

ss 



1.9 

51.4 
12.8 



7.4 
-2.3 
13.2 



5.0 
0.4 
5.6 



Table 6.11: A summary of the calculated mean uncertainties and systematic shifts in GeV for 
the M ± + M i endpoint analysis. FO denotes the results with fixed origin and NFO with 
non-fixed origin. 
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Figure 6.34: Theoretical kinematic endpoints for M ± + M t in GeV: The of coming from 

the f with the light quark and the p 1 coming from the X2 and the same light quark from 
Si,di — > s,d + X.2- 
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Figure 6.35: Example for the linear and the gaussian fit to the M ± + Mi distribution at 
mo = 181 and rriiM = 350. In this distribution 18714 events are represented. The triangle shows 
the theoretical kinematic endpoint at 768.8 GeV. 
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Figure 6.36: The difference of the theoretical and the measured M ± + M = endpoint in GeV 

Pi H H2 Q 

obtained with a linear fit for all 25 points. It shows a large shift of 51.4 GeV with a root-mean- 
square deviation of 12.8 GeV. 
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Figure 6.37: The theoretical M+„ + Mi„ endpoint minus the measured endpoint in GeV ob- 
tained with the linear fit after applying the constant correction of 51.4 GeV. The values decrease 
significantly for decreasing values of uiq. 
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Figure 6.38: The theoretical kinematic M+ n + -Mi endpoints as a function of the measured 
gaussian maximum for all 25 investigated points in the mo-Tti/2 plane. The gray fit with C = 1.9 
is forced to go through the origin whereas the black with C = 1.3 has an optimised ordinate 
value of 254.3 GeV. 
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Figure 6.39: The difference of the theoretical and the measured M ± + Mi endpoint in GeV 



obtained with the gaussian fixed origin method, 
systematic shift of 2.3 GeV 



The measured values are too high with a 
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Figure 6.40: The difference of the theoretical M + -\-M t endpoint and the measured endpoint 

r\H H2 * 

in GeV obtained with the gaussian non-fixed origin method. 
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6.2.10 The kinematic limit of M p ± p ^„ for bottom quarks 
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Indeed, with the b-tagging it is possible to select events with bottom quarks. But the top 
quark which is used in the sections 16.2.131 to 16.2.151 decays almost entirely into bottom quarks 
such that the b-tagging is also sensitive to top events. Nevertheless, there is a possibility to 
distinguish both kinds of events, as shown in figure 16.411 In the invariant mass distribution of 
both opposite-sign rho mesons versus the invariant mass of a single rho with the associated top 
or bottom quark, both kinds of events are well separated. Both processes are treated separately 
in the following. The analysis in this section is based on 25 data samples with between 4224 
and 4901 bottom quark events. 



For the kinematic limit of M p ± p ^i ) , which is shown in figure 16.421 for all 25 points, the 
gaussian method cannot be applied for two reasons. First, the shape of the distribution changes 
significantly within the investigated region. In some cases the region around the maximum of 
the distribution is flat, resulting in a large x 2 for the gaussian fit. Second, when applying the 
gaussian fit on the distribution the maximum value shows no linear dependence on the theoretical 
endpoint. This can be seen in figure IFTlTfl 



On the other hand, the linear fit (Figure I6.43|) provides a good possibility to measure the 
M p ± p ^ b endpoint. However, a large systematic correction (Figure I6.44J) is needed because there 
are about six to seven times more events with the lighter b\ than with the heavier 62 being 
responsible for the real kinematic endpoint. This can easily be understood by comparing tablc l6.il 
with table 16.41 The effect of this mixing is shown and explained in figure 16.47 



The same too 
dependence as in previous cases for the actually needed correction can be seen in figure 1^151 



Value 


Linear 


5statE P 1U 

s 

5S 


2.9 (2.9) 

50.3 (7.2) 

8.5 (8.6) 



Table 6.12: A summary of the calculated correction and mean uncertainties in GeV. The values 
in brackets correspond to the case of b\ being the input for the kinematic limit. 
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Figure 6.41: The invariant mass of both opposite-sign rho mesons versus the invariant mass of 
a single rho with the associated top or bottom quark in GeV. It is an example for the point 
I'. The left ellipse denotes the events with the bottom quarks whereas the top quark events are 
mainly located in the triangle on the right side. 




Figure 6.42: Theoretical kinematic endpoints in GeV of the invariant mass of both opposite-sign 
rhos, coming from the X2 and the f respectively, and the associated quark from 62 — > b + X.2- 
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Figure 6.43: Example for the linear and the gaussian fit to the M p ± p ^ b distribution at rriQ = 181 
and mi/2 = 350. The invariant mass distribution is based on 4498 bottom quark events. The two 
triangles show the theoretical endpoint values for events with b\ (532.9 GeV) and 62 (576.3 GeV), 
respectively. 
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Figure 6.44: The theoretical M p ± p ^ b endpoint minus the measured endpoint in GeV for all 25 
points. The mean value S is 50.3 GeV with a root-mean-square deviation of 8.5 GeV. If the 
b\ mass instead of the 62 mass is taken for the theoretical endpoint the value S changes to 
7.2 ± 8.6 GeV 
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Figure 6.45: The theoretical M p ± p ^ b endpoint minus the measured endpoint in GeV after ap- 
plying a constant correction of 50.3 GeV. The values decrease with decreasing values of itiq. 
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Figure 6.46: The theoretical kinematic M p ± p ^ b endpoints versus the measured gaussian max- 
imum for all 25 investigated points in the mo-mi/2 plane. Due to the flat shape near the 
maximum of one distribution, one point shows a large uncertainty. A unique mapping of the 
measured maximum to the theoretical endpoint can not be established since several theoretical 
endpoint values are associated with the same measured maximum. 
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Figure 6.47: The invariant mass distribution of both rho with the associated bottom quark. The 
bold black line shows 585 events coming from a 62, the bold gray line 3913 events coming form a 
61 and the light black line shows the superposition of both lines. Although the b\ events clearly 
dominate, the 62 events have a significant influence on the measurement of Ep M . The triangles 
show the theoretical endpoints at 533 GeV(6i) and 576 GeV(&2)- 



6.2.11 Two kinematic limits in M p ±„ for bottom quarks 
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The principles described in section 16, 2. Sl remain the same for this invariant mass distribution. 
In the following only the necessary changes and the results are described. The region for the 
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gaussian fit is changed from a symmetric area of 25% of the bin value which is used for the 
fitting centre to a symmetric area of 45% similar to section 16.2.81 The gaussian method for this 
distribution can be applied since the gaussian fit is stable for all 25 points and the assumed 
linear dependence is sensible (see figure Id".52|) . If the linear fit of Epf as a function of the 
fitted gaussian maximum Gi ± 5Gi is forced to go to the origin the resulting slope C is 2.3. The 
statistical uncertainty S s t a tEp on the endpoint measurement has then a mean value of 4.9 GeV. 
However, in this case a systematic shift S of —8.3 ± 14.1 GeV is needed in the whole region to 
improve the endpoint measurement (see figure l6*3H|) . If the fit is not restricted to intersect the 
origin the value C changes to 2.1, nearly no systematic shift S is needed and 5S decreases to 
6.0 GeV. The fit here intersects the ordinate at D = 21.6 GeV. The mean statistical uncertainty 
5 s tatEp M on the endpoint measurement decreases slightly to 4.3 GeV. The ratio between both 
slopes is 1.1. 

For the double linear fit method the second fit iteratively searches the area right of the 
second lepton peak where the slope is bigger than —0.25, instead of —0.2 like in section 16.2.81 
This value leads to a stable fit over the whole investigated parameter space. This method here 
again has the problem that a large shift of 46.9 ± 10.3 GeV is needed. However, it can be used 
to estimate the endpoint. Table l6~TTfl summarises the results. 



Value 



SstatEp 

S 

SS 



XT 



Linear Gaussian FO Gaussian NFO 



5.8 (5.8) 

46.9 (28.7) 

10.3 (7.7) 



4.9 
-8.3 
14.1 



4.3 
0.2 
6.0 



Table 6.13: A summary of the calculated mean uncertainties and systematic shifts in GeV for 
the second endpoint in M p ±t,. FO is with fixed origin and NFO with non-fixed origin. 



The endpoint for the first rho with the associated bottom quark is determined with a linear fit 
(Figure l6.56|) the same way as it is described in section lo.2.81 In figure lo351 the theoretical value 
is shown which corresponds to events coming from the 62. In figure l6~57l the large systematic 
shift S = 36.5 ± 3.7 GeV can be seen. There is no significant tuq and m\ / 2 dependence on that 
shift as shown in figure l6~58l This is similar to the first M p ± q endpoint for the light quarks 
(Figure EHHI)- 



Value 


Linear 


SstatE P M 
S 

8S 


2.6 (2.6) 

36.5 (-4.8) 

3.9 (3.7) 



Table 6.14: A summary of the calculated correction and mean uncertainties in GeV for the first 
endpoint in the M p ± q distribution. 
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Figure 6.48: Theoretical kinematic endpoints in GeV for the invariant mass of the second rho 
coming from the f and the associated bottom quark coming from the sbottom in the 62 
decay. The endpoints decrease significantly for decreasing values of uiq. 
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Figure 6.49: Example for the two linear fits and the gaussian fit to the M p ± b distribution at 
mo = 181 and min = 350. In this distribution 5636 events after the cut at 30% of the maximum 
invariant mass of both rhos (see section I6.2.8J) are represented. The two triangles show the 
theoretical endpoint values for events with 61 (244.7 GeV) and 62 (264.5 GeV), respectively. 
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Figure 6.50: The difference of the theoretical and the measured M ±, endpoint in GeV obtained 
with the linear fit for all 25 points. It shows a shift of 46.9 GeV with a root-mean-square 
deviation of 10.3 GeV. 
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Figure 6.51: The theoretical M ±. endpoint minus the measured endpoint in GeV obtained with 
the linear fit after applying the constant correction of 46.9 GeV. 
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Figure 6.52: The theoretical kinematic M n ± b endpoints as a function of the measured gaussian 



Pi 



maximum for all 25 investigated points in the mo-mi/2 plane. The gray line is forced to go 
through the origin whereas the black has an optimised ordinate value of 21.6 GeV. 
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Figure 6.53: The difference of the theoretical and the measured M +. endpoint in GeV obtained 

Pi ° 

with the gaussian fixed origin method. The measured values are too low with a systematic shift 
of -8.3 GeV 
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Figure 6.54: The difference of the theoretical and the measured M+, endpoint in GeV obtained 
with the gaussian non- fixed origin method. 
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Figure 6.55: Theoretical kinematic endpoints in GeV for the invariant mass of first rho coming 
from the X2 and the associated bottom quark coming from the sbottom in the 62 - ► b + X2 decay. 
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Figure 6.56: Example for the linear fit at m$ = 181 and m\/2 = 350 after applying a cut on the 
invariant mass M p ± b at the second endpoint. In this distribution 2036 events are represented. 
The two triangles show the theoretical endpoint values for events with b\ (514.7 GeV) and 62 
(556.4 GeV), respectively. 
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Figure 6.57: The difference of the theoretical and the measured M+, endpoint in GeV obtained 
with the linear fit. It shows a shift of 36.5 GeV with a root-mean-square deviation of 3.9 GeV. 
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Figure 6.58: The theoretical M n ± b endpoint minus the measured endpoint in GeV obtained with 



Pi- 



the linear fit after applying the constant correction of 36.5 GeV. 



6.2.12 The kinematic limit of M + + M=f„ for bottom quarks 
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For this kinematic endpoint the linear and the gaussian method can give a good estimate on 
its value. An example of a M_±. + M i, distribution is given in figure l6~HUl The data sample is 
the same as in section [6.2.101 

For the linear fit the large correction S = 89.3 ± 15.1 GeV is again caused by the origin of 
the bottom quark. A clear mo dependence can be seen in figure IO)21 which should be included 
into the correction. Similar to the light quark case the gaussian method can be applied as well. 
In order to achieve good results for the gaussian fit the symmetric area is chosen to be 35% 
of the maximum bin value. By comparing the two linear fits for the light quark gaussian case 
in figure 16.381 with the bottom quark gaussian case in figure 16.631 an interesting similarity can 
be found. In both cases the ratio between the slope in the fixed origin fit and the slope in the 
non- fixed origin fit has an universal value of 1.5. A summary of the measured values is given in 
table EH 
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Value 



SstatEp 
S 

ss 



IT 



Linear Gaussian FO Gaussian NFO 



3.2 (3.2) 
89.3 (37.2) 
15.1 (12.6) 



9.7 

6.2 

21.3 



6.5 
1.0 

6.4 



Table 6.15: A summary of the calculated mean uncertainties and systematic shifts in GeV for 
the M ± b + M = b distribution. FO is with fixed origin and NFO with non-fixed origin. 
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Figure 6.59: Theoretical kinematic endpoints for M ±. + M t, in GeV: The of coming from 



Pi 



Pi ' 



the f with the bottom quark and the p 1 coming from the X2 and the same bottom quark from 

the 62 — ► b + X2 decay. 
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Figure 6.60: Example of a linear and a gaussian fit to the M +.+M = h distribution at mo = 181 



and m 



Pi 



1/2 



350. In this distribution 4498 events are represented. The two triangles show the 



theoretical values for events with b\ (653.1 GeV) and 62 (706.1 GeV), respectively. 
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Figure 6.61: The difference of the theoretical and the measured M ±. + M t, endpoint in GeV 

Pi ° P2 

obtained with the linear fit for all 25 points. It shows a large mean shift of 89.3 GeV with a 
root-mean-square deviation of 15.1 GeV. 
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Figure 6.62: The theoretical M ± b + M i fe endpoint minus the measured endpoint in GeV ob- 



tained with the linear fit after applying the constant correction of i 
significantly for decreasing values of rriQ. 



.3 GeV. The values decrease 



800 




















700 














^ 


*-*" 


** 


600 














-JF 






£-500 
a> 
O 
"400 




















Q. 
LU 




















300 




















200 




















100 




, , i , , 




, , I , , 




, , I , , 


, , I , 


, , i , , 


, , I , 


°t 


)' ' 


50 


100 


150 


200 


250 


300 


350 


400 



G [GeV] 

Figure 6.63: The theoretical kinematic M +. + Mi, endpoints as a function of the measured 

P\ ° P2 

gaussian maximum for all 25 investigated points in the mQ-m\/2 plane. The gray line with 
C = 2.0 is forced to go through the origin whereas the black with C = 1.3 has an optimised 
ordinate value of 234.5 GeV. 



73 




Figure 6.64: The difference of the theoretical and the measured M ±, + M =, endpoint in GeV 

P\ P2 

obtained with the gaussian fixed origin method. The measured values are too low with a mean 
systematic shift of 6.2 GeV 
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Figure 6.65: The difference of the theoretical and the measured M ±, + M t, endpoint in GeV 

Pl ° H2 

obtained with the gaussian non-fixed origin method. 
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6.2.13 The kinematic limit of M p ± p ^ g for top quarks 
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For the analysis of events with top quarks it is necessary to use the endpoint formulae given 
in section 15.31 The input for the upper limit is the mass of the £2 and for the lower limit the 
lighter t\. For the following analysis 25 data samples with about 1750 events have been used. 
Similar to the previous cases, the gaussian method cannot be applied here which can be seen 
in figure l6~T01 The lower kinematic limit which is given by formula (|5.34|) does not provide any 
information which can be used for the mass reconstruction since the theoretical value hardly 
depends on m§ and mi/2- It is 193.5 ± 0.1 GeV over the whole investigated area. The upper 
limit given by formula (J5.33J) varies from 586.7 GeV at ttlq = 167 and m\/2 = 340 to 615.1 GeV 
at mo = 193 and m.1/2 = 360 which can be seen in figure l6~HHl The linear fit gives the endpoint 
which arises from the t\ mass with adequate precision. However, the real kinematic endpoint 
is given by the £2 mass. Thus, a quite large correction S of 155.5 ± 5.7 GeV is necessary. The 
reason for this is the stop mixing. The principle is described in section [(j.2.101 for the sbottom 
mixing. 



Value 


Linear 


5statE P M 

s 

5S 


3.1 (3.1) 

155.5 (0.4) 

5.7 (5.6) 



Table 6.16: A summary of the calculated correction and mean uncertainties in GeV. The values 
in brackets correspond to the case of £1 being the input for the upper kinematic limit. 
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Figure 6.66: Theoretical kinematic endpoints in GeV for the invariant mass of both opposite-sign 
rhos, coming from the X2 and the f respectively, and the associated top quark from £2 — ► t + Xi. 
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Figure 6.67: Example of a linear and a gaussian fit to the M p ± p ^ t distribution at rriQ = 181 
and mi/2 = 350. The invariant mass distribution is based on 1859 top quark events. The 
two triangles show the theoretical values for events with t\ (446.1 GeV) and £2 (601.2 GeV), 
respectively. 
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Figure 6.68: The theoretical M p ± p ^ t endpoint minus the measured endpoint in GeV. The mean 
value S is 155.5 GeV with a root-mean-square deviation of 5.7 GeV. If the t\ mass is taken for 
the theoretical endpoint the value S changes to 0.6 ± 5.6 GeV 
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Figure 6.69: The theoretical M p ± p ^ t endpoint minus the measured endpoint in GeV obtained 
with the linear fit after applying a constant correction of 155.5 GeV. 



77 



615 


I i 














i • — '< * ' 


610 




— • 


i • — 




i — • — i 








„605 

> 

0) 








— • — i 


O 


■ — • — i 


"600 

Q. 
UJ 

595 


'«* — \ 


i — • — 




■ — m^—r ' 
1 , , , 


,!,,,,!,,, 


590 


, , , I 




585 


i — • — i 


, I , 





240 



250 260 270 

G [GeV] 



280 290 



Figure 6.70: The theoretical kinematic M p ± p ^ t endpoints versus the measured gaussian maxi- 
mum for all 25 investigated points in the rriQ-mi/2 plane. A unique mapping of the measured 
maximum to the theoretical endpoint can not be established since several theoretical endpoint 
values are associated with the same measured maximum. 



6.2.14 Two kinematic limits in M p ±„ for top quarks 
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The theoretical lower limit of the first rho with associated top quark distribution is constant 
at a value 192.4 ± 0.8 GeV and for the second it is 176.3 ± 0.7 GeV. The upper limit of the 
second distribution shown in figure l6~TTl varies from 204.7 GeV at mo = 167 and mi/2 = 340 to 
376.0 GeV at mo = 193 and m\/2 = 360. The linear fit method can in principal be applied in 
the area with a too value around 170 GeV. However, for too values around 190 GeV this is not 
possible, even if the cut at 30% of the invariant mass of both rhos is used. An example for the 
first case at mo = 167 and m\/2 = 340 is given in figure l6~72*l and figure IB~73*1 shows an example 
of the latter at mo = 187 and m\/2 = 340. The gaussian method is applicable here as can be 
seen in figure l6~7H The only necessary change is that due to the large top mass the origin in 
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the Ep T -G plane used in the previous sections has to be moved from (0,0) to (175,175). This is 
called hence the "adapted" origin. The formula (|6.8j) has to be changed to 



Epf = C ■ (d - 175) + D + 175. 



(6.14) 



If the adapted origin is fixed the slope C is 5.6. If it is not, it decreases to C = 4.5 with D = 15.4. 



Value 


Gaussian FO Gaussian NFO 


SstatE P M 
S 

5S 


3.7 2.9 

-14.6 -1.2 

21.6 6.8 



Table 6.17: A summary of the calculated mean uncertainties and systematic shifts in GeV for 
M p ± t endpoint measurements. FO denotes the values with fixed adapted origin and NFO with 
non-fixed adapted origin. 



The upper limit of the first distribution shown in figure 16.771 varies from 555.0 GeV at 
mo = 187 and mi/ 2 = 340 to 610.1 GeV at tjiq = 173 and m 1 / 2 = 360. In order to achieve 
stable results the linear fit area has to start at the bin with 60% and ending at the bin with 5% 
of the maximum bin content. This change is necessary, because - due to the much lower statistics 
- for some points in the parameter space MINUIT [22] tries to fit a single bin. As can be seen 
in figure |p~T51 the linear fit actually measures the endpoint given by the t\. Therefore also the 
shift S and its uncertainty for measuring this endpoint is given in table 16.181 Figure 16.791 and 
figure IdTHOI however, show only the results for the real kinematic endpoint which is determined 
by i 2 . 



Value 


Linear 


5statE P M 

s 

6S 


3.2 (3.2) 

138.8 (-5.7) 

12.2 (4.2) 



Table 6.18: A summary of the calculated mean uncertainties in GeV for the p 1 endpoint in 
M p ± t . The values in brackets correspond to the case of t% being the input for the kinematic 
limit. 
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Figure 6.71: Theoretical kinematic endpoints in GeV for the invariant mass of the second rhos 
coming from the f and the associated top quark coming from the stop in the ii 
The endpoints decrease for decreasing values of ijiq. 



t + X.2 decay. 
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Figure 6.72: Example of M ± t at rriQ = 167 and m 



l l/2 



340 after the cut on M. 



±o=F 



P^P 



at 30% of 



its maximum invariant mass (see section 16. 2. 8|) . 2676 events are represented. The two triangles 

show the theoretical values for events with t\ (190.4 GeV) and £2 (204.7 GeV), respectively. In 

this case it is possible to measure the endpoint of the M+. distribution. 
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Figure 6.73: Example of M p ± t at m® = 187 and m-i/2 = 340 after the cut at 30% of the maximum 
invariant mass of both rhos. 2140 events are represented. The two triangles show the theoretical 
values for events with t\ (291.0 GeV) and £2 (371.8 GeV), respectively. In this case, even with 
the cut, it is not possible to measure the endpoint of the M ± , distribution. 



380 
















v*~ 


360 
340 


















320 
£300 

CD 

£280 
"£■260 






A, 


y 


/* 


Sk*"" 






240 


















220 
200 




3^ 














180 
160 


• 


I , , 


, , I , 




i i 




I , , , 


, I , , 



180 



190 



200 
G [GeV] 



210 



220 



Figure 6.74: The theoretical kinematic M ± t endpoints as a function of the measured gaus- 
sian maximum for all 25 investigated points in the rriQ-mi /% plane. The gray line is forced to 
go through the adapted origin whereas the black has an optimised adapted ordinate value of 
15.4 GeV. 
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Figure 6.75: The difference of the theoretical and the measured M n ± t endpoint in GeV obtained 



p 2 , 



with the gaussian fixed adapted origin method for all 25 points. The value is systematically too 
low on average by 14.6 GeV. 
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Figure 6.76: The difference of the theoretical and the measured M+. endpoint in GeV obtained 
with the gaussian non- fixed adapted origin method. 
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Figure 6.77: Theoretical kinematic endpoints in GeV for the invariant mass of first rho coming 
from the X2 and the associated top quark coming from the stop in the ti — > t + X2 decay. 
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Figure 6.78: Example for the linear fit to the M p ± t distribution at tuq = 181 and m 



1/2 



350 



after applying a cut on M p ± t at the theoretical second endpoint value. In this example 490 events 
are represented. The two triangles show the theoretical values for events with t\ (429.2 GeV) 
and £2 (582.3 GeV), respectively. 
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Figure 6.79: The difference of the theoretical and the measured M+. endpoint in GeV obtained 
with the linear fit for all 25 points. It shows a shift of 138.8 GeV with a root-mean-square 
deviation of 12.2 GeV. 
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Figure 6.80: The theoretical M+. endpoint minus the measured endpoint in GeV obtained with 
the linear fit after applying the constant correction of 138.8 GeV. 
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6.2.15 The kinematic limit of M t ± n + M =f„ for top quarks 
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Both methods can give an adequate estimate on the kinematic endpoint of the M+. + Mi, 

^^^_^ Pi r P2 

distribution. An example of such a distribution is shown in figure HTH21 The data sample is the 
same as used in section 16.2.131 The large correction S = 127.3 ± 28.5 GeV for the linear fit is 
mainly caused by the stop mixing. A significant thq dependence which should be included into 
the correction can be seen in figure HD£H Similar to the light and bottom quark case the gaussian 
method can be applied, which can be seen in figure IfTHol However, the gaussian function does 
not fit well in every case. In order to achieve adequate results for the gaussian fit the symmetric 
area is chosen to be 15% of the maximum bin value. The origin in the Ep T -G plane has to be 
moved from (0,0) to (350,350) since two top quarks are involved. This is hence again called the 
"adapted" origin. Thus, the formula (|6.8|) changes to 



Epf 



C -{Gi- 350) + D + 350. 



(6.15) 



The ratio between the slope of the fixed adapted origin fit and the slope in the non-fixed adapted 
origin fit has a value of 1.7. A summary of the measured values is given in table HQU1 which 
show that the ti definitely has an important influence on the endpoint measurement. 



Value 



SstatEp 

S 

ss 



w 



Linear Gaussian FO Gaussian NFO 



3.0 (3.0) 
127.3 (-52.6) 

28.5 (18.2) 



14.9 
1.3 

21.6 



1.9 
8.1 



Table 6.19: A summary of the calculated mean uncertainties and systematic shifts in GeV for 
endpoint measurements in M ± t + M i r FO is with fixed origin and NFO with non-fixed origin. 

The values in brackets correspond to the case of t\ being the input for the kinematic limit. 
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Figure 6.81: Theoretical kinematic endpoints for M ± t + M = t in GeV: The p^ coming from 



f'T 



-,o 



Pi 



the f with the top quark and the p 1 coming from the X2 and the same top quark from the 
£2 — ► t + X2 decay. 
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Figure 6.82: Example of a linear and a gaussian fit to M +.+M T , at mo = 181 and mi/2 = 350. 

Pi ^ P2 ' 

In this distribution 1859 events are represented. The two triangles show the theoretical values 
for events with t\ (578.6 GeV) and £2 (761.4 GeV), respectively. 
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Figure 6.83: The difference of the theoretical and the measured M ±. + M^. endpoint in GeV 
obtained with the linear fit. It shows a large mean shift of 127.3 GeV with a root-mean-square 
deviation of 28.5 GeV. 
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Figure 6.84: The theoretical endpoint minus the measured M+. + M r ^, endpoint in GeV ob- 
tained with the linear fit after applying the constant correction of 127.3 GeV. The values decrease 
significantly for decreasing values of mo- 
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Figure 6.85: The theoretical kinematic endpoints as a function of the measured gaussian maxi- 
mum for all 25 investigated points in the mQ-m\j2 plane. The gray line with C = 4.3 is forced 
to go through the adapted origin whereas the black line with C = 2.4 has an optimised adapted 
ordinate value of 169.0 GeV associated to the real ordinate value of —335.0 GeV. 
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Figure 6.86: The difference of the theoretical and the measured M +, + Mi, endpoint in GeV 
obtained with the gaussian fixed origin method. 
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Figure 6.87: The difference of the theoretical and the measured M+. + Mi. endpoint in GeV 
obtained with the gaussian non-fixed adapted origin method. 
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6.2.16 Examples of measurements with 30 fb 



-i 



In the following a Monte Carlo sample at I' is used which corresponds to an integrated lumi- 
nosity of 30 fb^ 1 if all efficiencies and acceptances were 100%. The aim is to test the methods 
introduced in the previous sections for a low statistics sample. Table 16.201 shows the results of 
the endpoint measurements and figure loTHSl to l6~100l show all associated distributions. 



Distribution 


Linear 


Gaussian FO 


Gaussian NFO 


Ep 1 


Mp P 


83.t> ±L8 


81.0*- 


81-2^ 


85.8 


M ppq 
M P2q 
M Pig 
EM M 


614.8^ 

320.0+1 3 7 
594.3g;l 

742.7g 2 7 8 


n.a 
300.2%H 

n.a. 
77D fi^ 15 - 3 


n.a. 
296.7± 2 7 ; 7 

n.a. 
77 A i± 10 - 3 


627.8 
288.0 
606.0 

768.8 


Mp pb 
M P2b 
M Plb 
J2M Pxb 


560.4|^' 

91Q 9 ±25.4 

542.3^ 9 4 
700.0±£ 7 


n.a. 
265.2|£° 

n.a. 
679.(^11:1 


n.a. 
262.4^° 6 

n.a 
686.3^ 


576.3 
264.5 

556.4 
706.1 


Mpp t 
M P2t 
M Pl t 


592.2^^ 

n.a. 
559.9±£° 

917.5±1 6 5 3 


n.a. 
285.9 ±216 

n.a. 
744 v±6l.2 


n.a 
290.2±«:1 

n.a. 
743.6I 34 ! 6 


601.2 
316.7 
582.3 
761.4 



Table 6.20: The measured values Ep M % s ^ t for 30 fb~ 
denoted with "n.a.". All values are in GeV. 



If endpoints are not available they are 



With the 13 measured endpoints it is now possible to reconstruct the masses of the involved 
sparticles. For the mass reconstruction the minimisation package MINUIT within the ROOT 
framework [22] is chosen. The input to the minimisation are the measured endpoint values 
with their associated formulae from section 15.21 and 15.31 and the measurement uncertainties 

5 = ^J ' stat 2 + syst 2 . 

First, the measurements of the dilepton mass and the masses with the light quark with the 
smallest 5 have been chosen. Thus, five endpoints are used to measure four masses. The result 
(with the masses calculated by ISASUGRA in brackets) is: 



M 



d/si 

M fl 



787.5 ± 35.5 GeV (782.6 GeV) 

261.6 ± 25.5 GeV (265.5 GeV) 
157.4 ± 26.5 GeV (150.1 GeV) 
144.3 ± 25.7 GeV (138.1 GeV) 



(6.16) 



The uncertainties given in 16.161 are noticeably larger if only the first four endpoints are used. 
There are two reasons for these, nevertheless, large uncertainties. First, the endpoints cannot 
be measured as precisely as it is possible for low tan f3 models due to the presence of neutri- 
nos in the decay chains. Second, the equations for the endpoints are strongly correlated and 
therefore do not contain sufficiently complementary information which could compensate the 
endpoint uncertainties. By adding the endpoint measurements which involve the bottom quark 
the uncertainties can be improved and the sbottom mass can be measured. Here nine measured 
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values are taken to determine five unknown masses. 



M 



d/si 
M i 2 

M x% 

M fl 

Xi 



778.4 ± 28.2 GeV (782.6 GeV) 

710.8 ± 25.7 GeV (724.9 GeV) 

256.9 ± 20.4 GeV (265.5 GeV) 
150.9 ± 21.2 GeV (150.1 GeV) 
138.2 ± 20.5 GeV (138.1 GeV) 



(6.17) 



With adding the endpoint measurements which involve the top quark the uncertainties on the 
reconstructed masses slightly increase due to the large endpoint uncertainties. However, the 
stop mass can be estimated. Here 13 measured values are taken for the determination of six 
unknown masses. 



M 



d/s\ 

M h 

M x% 
M fl 

X\ 



795.2 ± 29.9 GeV (782.6 GeV) 

743.1 ± 28.8 GeV (747.9 GeV) 

727.2 ± 27.2 GeV (724.9 GeV) 

272.4 ± 21.7 GeV (265.5 GeV) 

164.5 ± 22.6 GeV (150.1 GeV) 
151.9 ± 21.9 GeV (138.1 GeV) 



(6.18) 



The uncertainties on the masses are larger than the mass differences within the investigated area 
of the mSUGRA parameter space. However, it is at this level possible to give good estimates 
on these six sparticle masses. 





Figure 6.88: The invariant mass distribution 
of both opposite-sign rhos is based on 3062 
events. The triangle shows the theoretical 
endpoint value at 85.8 GeV. 



Figure 6.89: The invariant mass distribution 
of both rhos and the associated light quark is 
based on 2282 events. The triangle shows the 
theoretical kinematic endpoint of M p ± p ^ q at 
627.8 GeV. 
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Figure 6.90: The invariant mass distribution 
of the second rho and the associated light 
quark with the first rho background is after 
the cut (55%) based on 844 events. The tri- 
angle shows the theoretical endpoint value at 
288.0 GeV. 



Figure 6.91: The invariant mass distribution 
of the first rho and the associated light quark 
is after the cut on the invariant mass M p ± q 
at the second endpoint based on 1199 events. 
The triangle shows the theoretical endpoint 
value at 606.0 GeV. 
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Figure 6.92: Example for the linear and the 

gaussian fit to the M + + M =c distribution 
P\ i P2 y 

at rriQ = 181 and m 1 / 2 = 350. In this dis- 
tribution 2282 events are represented. The 
triangle shows the theoretical kinematic end- 
point at 768.8 GeV. 



Figure 6.93: The invariant mass distribu- 
tion of both rhos and the associated bot- 
tom quark is based on 550 events. The two 
triangles show the theoretical endpoint val- 
ues for events with b\ (532.9 GeV) and 62 
(576.3 GeV), respectively. 
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Figure 6.94: The invariant mass distribution 
of the second rho and the associated bottom 
quark with the first rho background is after 
the cut (30%) based on 642 events. The two 
triangles show the theoretical endpoint val- 
ues for events with b\ (244.7 GeV) and b 2 
(264.5 GeV), respectively. 



Figure 6.95: The invariant mass distribu- 
tion of the first rho and the associated bot- 
tom quark is after the cut on the invari- 
ant mass M p ± b at the second endpoint based 
on 249 events. The two triangles show the 
theoretical endpoint values for events with 
61 (514.7 GeV) and b 2 (556.4 GeV), respec- 
tively. 
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Figure 6.96: Example of a linear and a gaus- 
sian fit to the M ±, + M 1 , distribution at 

p 1 b pjb 

mo = 181 and m\i 2 = 350. In this distribu- 
tion 550 events are represented. The two tri- 
angles show the theoretical values for events 
with 61 (653.1 GeV) and b 2 (706.1 GeV), re- 
spectively. 



Figure 6.97: The invariant mass distribution 
of both rhos and the associated top quark 
is based on 230 events. The two triangles 
show the theoretical values for events with t\ 
(446.1 GeV) and i 2 (601.2 GeV), respectively. 
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Figure 6.98: The invariant mass distribu- 
tion of the second rho and the associated 
top quark with the first rho background is 
after the cut (30%) based on 280 events. 
The two triangles show the theoretical val- 
ues for events with i\ (257.1 GeV) and £2 
(316.7 GeV), respectively. 



Figure 6.99: The invariant mass distribution 
of the first rho and the associated top quark 
is after the cut on M p ± t at the theoretical 
second endpoint value based on 112 events. 
The two triangles show the theoretical val- 
ues for events with t\ (429.2 GeV) and £2 
(582.3 GeV), respectively. 
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Figure 6.100: Example of a linear and a 

181 
350. In this distribution 230 



gaussian fit to M+, + Mi, at mo 



and mi/2 

events are represented. The two triangles 
show the theoretical values for events with 
h (578.6 GeV) and i 2 (761.4 GeV), respec- 
tively. Especially this figure shows the ne- 
cessity to further optimise the fitting meth- 
ods for low statistic samples. 
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Chapter 7 

Conclusions 



If nature is described by mSUGRA in the investigated region around I' there will be some 
difficulties to reconstruct the investigated sparticle masses in the dileptonic decay of the X2. 
However, with the Large Hadron Collider the CMS and the ATLAS experiments would be able 
to find signals for Supersymmetry and to give an estimate of the sparticle masses. The results 
of this study can be summarised as follows: 

1 . Since the cascade decay of the X2 is dominated by tau leptons at high tan (3 the dileptonic 
invariant mass distribution has no sharp edge at the kinematic limit. 

2. The t — ► p v T channel offers a good opportunity to search for endpoints. 

3. In the investigated parameter space it is possible - at the Monte Carlo level - to reconstruct 
five to six sparticle masses in a 30 fb^ 1 data sample, with an uncertainty of about 20 to 
30 GeV, by using the methods introduced in this study. 

4. There is a linear dependence between the gaussian maximum and the theoretical endpoint 
for several mass distributions. Thus, endpoint measurements with a linear fit at the tail 
of the distribution can be well improved. 

5. In principle it is possible to obtain two kinematic limits - the limit of M ± and M ± - by 
measuring the M p ± q distribution with a cut on the invariant mass of both opposite-sign 
rho mesons. 

6. Top and bottom quark events can be well separated in the M pp — M pq plane. 

7. In the investigated parameter space the endpoints of the distributions involving bottom and 
top quarks are determined by the 62 and £2, respectively. However, the measured endpoints 
with the linear fit method correspond more to the b\ mass and t\ mass, respectively, 
because b\ and t\ significantly dominate the distributions. Thus, it is generally very 
important to take the mixing of the involved particles into account. 
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Chapter 8 

Outlook 



In this study many simplifications have been applied. In order to obtain a more realistic pre- 
diction, at least the improvements given below have to be included. Furthermore, even at the 
Monte Carlo level additional studies like the investigation of the 3-prong channel may give useful 
information. The next steps towards a more detailed analysis are: 

1. Estimate of background 

This study does not take Standard Model backgrounds into account. Additionally, there is 
also a mSUGRA background coming from taus in chargino and other decays. An estimate 
of these backgrounds is necessary. 

2. Estimate of detector effects 

Since this study is purely at the Monte Carlo level, no detector effects have been taken into 
account. The CMS programs OSCAR for the full simulation, ORCA for the reconstruction 
and FAMOS for a combined fast detector simulation and reconstruction are necessary to 
give an estimate of detector effects. 

3. Investigation of the 3-prong channel 

Not only the r — ► p v T channel can be used for the X2 discovery at high tan/3. Due 
to the nearly similar particles involved in the decay chain, the r — * it tt^tt v t channel 
should have related properties. Assuming the same efficiency and acceptance, its use would 
increase the statistics by a factor 2.5. 

4. Possible effect through spin correlations 

In the simulations no spin correlations are taken into account. However, it is possible to 
make an estimate how the spin correlations might have an effect on the invariant mass 
distributions. Since the statistics in the region of an kinematic endpoint depends on 
whether that special configuration is favoured or not the investigation of that configuration 
with respect to the spin gives a hint of the a possible effect. 
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Appendix A 

Problems using PYTHIA 6.220 



A.l Bug in the hadronic tau decay 

We found that in PYTHIA 6.220 the invariant mass distribution of the visible particles in the 
t^ — ► p°7r ± ^ r decay has a r-energy dependence which is unphysical. Further investigations 
showed that the energy distribution of the v T in the rest frame of the t^ depends on the r- 
energy in the laboratory frame of the program, which is a clear violation of Lorentz Invariance. 
According to the author the intended neutrino spectrum correction had not been obtained for 
taus with E T > 20m T , which is the limit where taus are boosted to their rest frame for better 
numerical precision [23]. The bug has been fixed in version 6.222, which has been released on 
the web [23] on 21 st January 2004. However, the author remarks that the correction of the 
neutrino spectrum is still crude. 

A. 2 Particle listing 

In the particle listings of PYTHIA 6.220 interfaced with ISASUGRA 7.69 we found some prob- 
lems depending on the value of MSTP(128). It is not clear whether the problem is caused by 
PYTHIA in this case. In the following the change of the particle listing due to different values 
of MSTP(128) is briefly presented. Particles in the documentation section have status three and 
particles in the rest of the particle listing have status one or two. 

1. MSTP(128) = 0: For this value there is a huge overlap between the documentation section 
and the rest of the particle listing. All intermediate states are declared here, but it happens 
that particles point to their children, but the children point to the mother particle of their 
mother. This causes huge problems when selecting special decay chains. 

2. MSTP(128) = 1: This is the value which is chosen in this study. Even here there are 
difficulties in finding the right decay chain, but it is possible to make unique assignments. 

3. MSTP(128) = 2: In this case the problem is that particles in the particle listing with status 
one or two do not always point to the right mother particle in the documentation section 
with status three. It happens that intermediate states are totally missing sometimes. Since 
the particles in the documentation section carry no information about their children, it is 
therefore difficult to distinguish between different decay chains. 

A.3 PYTHIA 6.220 and TAUOLA 

We tried to use TAUOLA to solve the tau problem described in section IA.1I but we did not 
achieve to run TAUOLA with PYTHIA and ISASUGRA. For Standard Model processes we 
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found no problem. C. Biscarat mentions this problem on her homepage [24] and we assume that 
the mixed order of particles in PYTHIA (see lA.2)) is responsible for that. 

A. 4 Calculated cross sections 

It seems that there is a dependence of the cross section values on the selected production 
channels. We found changes in the production cross section of about 20 %. 
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